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Abstract 


Ejection and prapagation of the particles which excite abnormal ionization in 
the polar upper atmosphere are investigated from some statistical considerations 
on ionospheric data f/-min and cosmic radio observation by riometer. 

Relations between the flares and the polar-cap blackouts insist the trapping 
of particles into the cloud which would be responsible for the geomagnetic dis- 
turbances; that is, a part of the energetic particles ejecting from the flare 
propagate with relativistic velocity in interplanetary space along the twisted 
solar magnetic line of force extending to the earth and excite the S-type polar- 
cap blackouts having sudden onset within a few hours after flare; on the other 
hand, another part of the energetic particles ejecting from the flare is trapped 
into the cloud and is carried with rather low velocity, then the particles would 
reach the earth by leaking from the cloud and excite the G-type polar-cap black- 
outs having gradual onset about 10 hours before commencement of geomagnetic 


storms. 


1. Introduction. ~~ 


Recently, new evidence of the high energy solar particles has been found from 
several studies. Anderson, et al. (1959) showed difinitely the existence of solar protons 
by means of balloon observations at Churchill. On the other hand, Shapley, et al. 
(1958) and Haura, et al. (1958) (1959) showed from analysis of f-min ionosphere data 
that the intense polar blackouts caused by enhanced inoization in lower ionosphere at 
polar cap regions occurs after an intense solar flare associated with type IV solar radio 
outbursts. , 

These evidences have also been identified by riometric observation of the galactic 
radio noise (Little et al. 1959, Reid et al 1959, Hultqvist 1959) and by observation of 
VHF polar circuits (Bailey 1957, 1959). 

Many papers on this subject have now been presented from the experimental and 
theoretical points of view. Lately, Obayashi and Hakura (1960) pointed out the tend- 
ency that the polar blackout particles originating to the west of solar central meridian 
reach to the earth earlier than those from the east, and they explain this tendency 
under the conception that the particles tend to travel easily and speedily along the 
curved solar magnetic lines of force which are twisted towards the east (Biermann 


1957, Parker 1958). 
€ i*) 
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Here, the ejection and propagation of the solar particles responsible for the polar 
blackout events has been discussed upon the consideration of both the corpuscular 
cloud which moves along the straight path from the sun at a low speed and the parti- 
cles which moves along the twisted solar magnetie lines of force at a high speed. 


2. Sudden and Gradual Onset of Polar-Cap Blackouts 


During the period from January, 1956, through December, 1960, fifty pre-sse polar- 
cap blackouts were observed by means of f-min ionosphere data at Thule, Greenland 
(Geograph. 77.29N, 69.10W; Geomag. 89.0, 357.8, near the geomagnetic north pole) 
and at Resolute Bay, Canada (Geograph. 74.41 N, 94.55 W ; Geomag. 88.9, 289.3, near 
the magnetic north pole), and also by riometric observation at various stations (Reid 
et al. 1959). 


Sudden Polar Absorption 


Sth Jul.19s7 


Flare (11+) 


ith Feb. 1958 


Fig. 1 Typical examples of polar blackout events observed by f-min of 

ionospheric data at Thule. The upper figure is Sudden-type event which 

started on 3rd July, 1957 and the bottom figure is Gradual-type event 
which started on 10th February, 1958. 


The relevant details of these events have been compiled in Table I with their 
associated solar and magnetic events. 

After looking over the development of the polar blackouts, we easily come to the 
conception that there are two types of polar-cap ionization events, one having a sudden 
onset within two hours (S-Type) and the other having a gradual onset (G-Type) (Sinno 
1916 a). 


The examples of S- and G-types of the polar blackouts observed from Ff-min ionos- 
pheric data at Thule are shown in Fig. 1. 


3. Propagations of Polar Blackout Particles 


3. 1 Propagation times between polar blackouts and geomagnetic disturbances 

As has been investigated (Hakura et al. 1959), we can see from Table I that almost 
all cases of the polar blackout events are preceded, the same as the magnetic storms 
(Sinno et al. 1958), by type IV solar radio outbursts. Up to eighty percent of polar 
blackouts correspond to type IV outbursts from the table, The time delays of the 
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pre-ssc polar blackouts (4t) and of the sudden commencements of geomagnetic storms 
(AT) after the flare are also tabulated in this table. 


8 8 €88ssz 
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At 


St (Hours) 


Heliograph itude 
iegraphic Longi AT (Hours) 


Fig. 2 Relation of the propagation time Fig. 3 Heliographic longitudinal distribution 
between polar blackouts st and magnetic of the propagation times of polar blackout 
storms 47. Polar-cap blackouts with events. 


sudden onset and gradual onset events 
are indicated by © and @, resrectively, 
and uncertain events indicated by O. 


The relation beween 4¢ and 4T is shown in Fig. 2, indicating S.and G-Types 
events as double circles and black circles, respectively, and uncertain type events as 
white circles. Remarkable characteristics we could easily deduce from Fig. 2 are the 
existence of two wings; one wing contains almost all S-type events starting within a 
few hours in 4¢ and the other wing contains almost all G-type events starting about 
ten hours before ssc (Sinno 1961a). 

These characteristics lead to the conception that the latter shows some inter-rela- 
tion with the corpuscular clouds responsible for the geomagnetic distubances. It may 
be suggested that a part of the energetic solar particles emitted from the flare is 
trapped into the cloud having a intense magnetic field and carried with it and then 
impinges into polar atmosphere by gradual leakage from the cloud. 


3. 2 Heliographic longitudes 

The time delays of the polar blackouts after the solar flares 4t are plotted in Fig. 3 
according to their heliographic latitude. As was pointed out (Obayashi et al. 1960), 
the polar-cap blackouts are excited by the flares on the eastern side. These observed 
time-delays and western excess of the polar-cap blackouts are favorable to the concep- 
tion that the solar magnetic line of force are twisted towards the east owing to the 


4 K. SINNO 


solar rotation, since the particles easily propagate along the line of force. However, 
there are several events which belong almost to G-type as shown in Fig. 3, the trav- 
elling time corresponding to many hours, so we can not simply consider that the 
particle travels along the line of force from the sun, but it may travel ina complicated 
manner, e. g., the particle is trapped within the corpuscular cloud moving slowly from 


the sun before it leaks from the cloud. 


Number of Occurrence 
Absorption in DB 


Heliographic Longitude Heliographic Longitude 


Fig. 4 Heliographic longitudinal occur- Fig. 5 Maximum absorptions of polar black- 
rences of the polar blackouts. outs (estimated from riometric observations) 


vs heliographic longitude. 


Longitudinal occurrences of the flares responsible for polar blackouts are shown 
in Fig. 4. We can find the significant west hemisphere excess in respect of S-type 
and also solar cosmic ray increase, as has been presented by Reid and Leinbach (1959) 
etc., but we see G-type blackouts correspond to the flare of rather wide spread longi- 
tudes. 

We know that trapping of the energetic particles into the cloud is more effective 
by the flare of eastern heliographic longitude than those of western heliographic longi- 


Number ot Blackout 


<———_ Time of Flare ( Average) 


Storm Time (Hours) 
Fig. 6 Number of occurrences of polar blackouts with respect to the 
geomagnetic storm time. 
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tude, as shown by the present author from investigation of cosmic ray storms (Sinno 
1961b). 

If G-type blackout is evidence of leakage of the trapped energetic particles in the 
cloud, the widespread heliographic occurrence of G-type would be caused by combina- 
tion of both trapping and leaking tendencies, predominating to east and west hemi- 
sphere, respectively. 

Here, we would like to point out the fact that ionosphere absorption of galactic 
radio noise during polar blackout does not show the heliographic western excess, as 
show in Fig. 5. It is noticeable to know that G-type polar blackout is usually more 
intense than S-type, so that energy spectrum of G-type is much steep comparison to 
that of S-type. Therefor, it is plausible to consider that the particles of G-type black- 
out had been trapped into the cloud. 


3. 3 Storm time 

' Number of occurrence frequency of polar blackouts with respect to the geomagne- 
tic storm time takes the maximum at the time of near ssc, as shown in Fig. 6. 
Though, this tendency is an average of all polar blackouts, it may represent only a 
tendency of G-type, since G-type polar blackout is usually more intense than S-type 
in low energy range. The fact may show that the lower energy particles of solar 
cosmic ray responsible for polar blackouts are easily trapped into the cloud, compared 
with the higher energy particles of solar cosmic ray which reach their mamimum 
within a few hours after the flare, because the trapping ability may be inversely pro- 
portional to energy of the particles as has been discussed by the present author (Sinno 
1961b). Thus, S-type polar blackout with sudden onset (~100 MeV as proton) might 
be a trace of unusual increase of solar cosmic rays since. they have rather flat spectrum 
towards the higher energy. Actually, we have always observed S-type polar blackouts 
at the time of unusual solar cosmic ray increase. 


3.4 Ejections of Polar Blackout Particles 

The fact that the polar-cap blackout phenomenon iis excited by the major solar 
flare associated with the type IV radio emission was already confirmed statistically 
on the IGY data (Hakura et al. 1959, Thompson et al. 1960). It is also well known 
that the cause of the type IV solar radio outburst is synchrotron radiation from the 
relativistic electrons in the solar outer corona (Boischot et al. 1957). The existance 
of a large number of those energetic electrons leads us to the conception that there 
are also many protons with relativistic energy. 

When these protons are ejected from the sun, the corpuscular cloud responsible 
for the geomagnetic storm’ and also the cosmic-ray storm should be ejected, since it 
is true that both the storms are excited by the solar flare associated with the type IV 
radio noise. 

As we already known, big solar radio outbursts are usually consisted by two parts, 
pre-flare (first part) and post-flare (second part) outburst (Dodson et al. 1953). Pre-flare 
outburst.is characterized by abrupt onset of outburst having very wide frequency ranges 
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and closely associated with SJD phenomenon of ionospheric radio propagation. On 
the other hand, post-flare outburst called type IV, is characterized by gradual onset of 
outburst having rather low frequency below several hundred Mc, and closely connected 
with the cloud responsible for the geomagnetic storm and the cosmic ray storm (Sinno 
et al. 1958). 

In order to answer the question of which part of the outbursts is evidence of the 
polar blackout particles, we will point out the following observational fact. We have 
several events where comic ray increase and S-type of polar blackout started before 
onset of post-flare outburst. That is, 23rd February, 1956, pre-flare outburst started 
at 03:35, cosmic ray increase started before 03:50 and post-flare outburst started after 
04:00. And, 12th November, 1960, pre-flare outburst started at 13:21, cosmic ray 
increase and polar blackout started at about 13:40 and post-flare outburst started 
after 14:00. 

On the other hand, we have also several events where G-type of polar blackout 
was excited appreciable pre-flare outburst (11th September. 1957 and 9th February, 
1958). 

Furthermore, we have a statical investigation which supports this correspondence, 
that is, the polar blackout having short propagation time corresponds to big solar 
outburst on higher frequencies (Sakurai et al. 1961), and to solar outburst with simul- 
taneous onset over wide band (Sinno et al. 1958). The outburst having these charac- 
teristics should be the pre-flare outburst itself. 

Then, we come to the conclusion that pre-flare outburst is the evidence of ejection 
of solar energetic particles and post-flare outburst is the evidence of trapping of the 
particles, though position of the flare may also have an important effect on excitation 
of polar blackout. 


Table II Courses of Polar Blackouts 


Polar Blackout 
(Phenomenon) 


Solar Outburst Interplanetary Space 


Flare (Location) (Ejection) (Fropagation) 


M Along to Solar 
West . Pre-Flare Magnetic Field S-Type 
Hemisphere : 


Post-Flare Trapped into Cloud G-Type 


East 
Hemisphere 


. Along to Solar 
Pre-Flare Magnetic Field No Blackout 


See ee ee esta TO 

Courses of development of polar blackouts from solar flares are shown in Table II. 
If pre-flare solar radio outburst occurred ‘in western heliographic longitudes, ejected 
energetic particles propagate along to the twisted solar magnetic line of force and 
could excite S-type of polar blackout and cosmic ray increase with short time delays, 
however if the same pre-flare outburst occurred in eastern or central heliographic 
longitudes, ejected particles propagate also along to the field and could not strike to 
the earth, then no polar blackout is excited. On the other hand, G-type polar blackout 
is excited by post-flare outburst occurring anywhere in visible heliographic longitudes, 
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since lateral dimension of the cloud is quite large as several A. U. near the earth orbit. 


First part 


| 

| 
Solar Radio < Second part 
Qut burst o= 


= 
! 
j 
| 
I unusual increase 
| 
l 
| 


Cosmic Ray 
Storm p Forbush decrease 


Polar Cap 
Blackout 


Geomagnetic 
StormH. 


Fig. 7 Schematic representations of the two families of the earth storms. 


Concluding Remarks 


Ejection and propagation mechanisms of polar blackout particles are investigated 
from some statistical considerations. 

It may be said that there are two ways of ejection of the energetic particles through 
association with either pre-flare or post-flare solar radio outbursts. The former particles 
having energy up to cosmic ray energies propagate speedily (~10'°cm/sec) along to 
solar magnetic line of force twisted to the eastwards, and excite both S-type polar 
blackout with sudden onset and solar soft cosmic ray increase, if the earth comes 
across the path of the particles. On the other hand, the particles having comparatively 
low energies are trapped easily into the cloud and carried with rather lower velocity 
(~108cm/sec) in interplanetary space and particles leaking from the cloud excite G-type 
polar blackout with gradual onset, if the earth comes across the path of the cloud. 

Auroral particles during magnetic storm might be also a trace of those of lower 
energy ranges and impinge into the upper atmosphere only during a restricted period, 
because they are trapped strongly into the cloud by their low momentum. 

In concluding, we come to a conclusion that earth storms might be divided into 
two families, one constructed by solar cosmic ray including S-type of polar blackout, 
SFE of geomagnetism and SID of radio propagation, associated with pre-flare (first 
part) solar radio outburst, and the other constructed by cosmic ray storm (Forbush 
decrease), G-type of polar blackout including aurora, geomagnetic storm and ionospheric 
storm, associated with post-flare (second part) solar radio outburst, as shown in Fig. 7. 
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The author wish to express his thanks to the members of the Cosmic-ray Research 


Committee in Japan for their valuable discussions. 


References 


Anderson, K. A., Arnory, R., Peterson, L. and Winckler, J. R., (1959) J. Geophys. Res., 64, 1133. 


Bailey, D. K., (1957) J. Geophys. Res., 62, 431. 

Bailey, D. K., (1959) Proc. IRE, 47, 255. 

Biermann, L., (1957) Observatory, 77, 109. 

Boischot, A. and Denisse, J. F., (1957) Comptes Rendus, 245, 2194. 
Dodson, H. W. Hedeman, E.R. and Owren, L., (1953) Ap. J., 118, 169. 
Hakura, Y., Takenoshita, Y. and Otsuki, T., (1958) Rep. Ionos. Res. Japan, 12, 459. 
Hakura, Y. and Goh, T., (1959) J. Radio Res. Lab. Japan, 6, 635. 
Hultqvist, B., (1959) Tellus, 11, 332. 

Little, C.G. and Leinbach, H., (1959) Proc. IRE, 47, 315. 

Obayashi, T. and Hakura, Y., (1961) J. Radio Res. Lab. Japan, 7, 379. 
Parker, E. N., (1958) Ap. J. 128, 664. 

Reid, G.C. and Leinbach, H., (1959) J. Geophys. Res., 64, 1801. 
Sakurai, K. and Maeda, H., (1961) J. Geophys. Res., 66, 1966. 
Shapley, A. H. and Knecht, R. W., (1958) NBS Rep., No. 5596. 

Sinno, K., (1961a) J. Radio Res. Lab. Japan, 8, 17. 

Sinno, K. and Hakura, Y., (1958) Rep. Ionos. Res. Japan, 12, 285. 
Sinno, K., (1961b) J. Radio Res. Lab. Japan, 8, 289. 

Thompson, A. P. and Maxwell, A., (1960) Nature, 185, 89. 


Characteristics of Solar Energetic Particles which Excite Polar-Cap Blackouts 


(1S) 


euawoueyd pazeloosse YIM s}U9AI UOTZezIUO! deo sejod jo ysr] 


ar) AY | OR ie) TL | OEM Cen | Zee ie oe) = aa = a au((9)| 7 | st: 00 so 5 Ut SF ='00 20 
aa = 3 QD | = O€ | 0S: TEL OL| 34] XO OF | gees Te GE |e = 
8V | FS = |@-0l| +E | OH StS | OG: 60 2 | Xa | Zt 9 | 08: <¢ $6 | HE? © 48 | S¥:60 SZ ot | On: SE Se 
| ST | & | O8M ES | Sb= Siri) = = te a Hey § LG WeSOr BE | ri 
ahaigae re tt tt ee ae | aca va . $i - = BiG 2 66 | 00:01 IT SUL | feeb VE 
OL | 8@ AY | G0- te) +c | PEM StS | Sosa) OD |) ci — | 0: toh) Al S&S T€ | 0€: 20 OT s Gé - 10 IT 
= Iv AT |) 6060 tore: |. ere -OUSs | See Some — a aa ad a = aa Su | So: 10 9¢ 
vI | O€ AL | 9¥:9T| +€ | SEM 9S | H:9T 02 | 00) S Er 0020) tee) Lal cS oem | SS) Sy | ge ee Tv: Ge 12 
Pr. Les AL | 96-61 | .€ | S08 8S | 0°60 9210) | < GG. ipSl= Sa. 968) h— ia. = = s 91:00 6¢ 
9 | We AT | O€-€T| € | SOM EIN | 2E:E1 12 | 09) & TS 062 GE 1G Ww .L. |) “S ce | .St sleet s SV: &I 2 
0€ | OF AI | 00:0) € | SOA LIN | Sh:ZO IL | OD |s°0 so |) 00. CL clits 2 6I | 0€: 60 @L s 94:00 €I 
8 87 AI | OT:€E| € | 9EM GZS | €L:e1 201 09]| 6 c& | 00:12 co | — —_ | Penuyuos s 00: €I ¥0 
hiook 8€ AI | 60: €1 € | SOM 9EN | 6S: cI IE | — a - io Hw} S OOT | ST: 71 TE s VI: £0 20 
8°61 | VE AL | O&- GO) +e | S6a CES | 7:60 82)) OO | 6 £6 | 00: €i6e ) ML} SS Aer SO 66 Su | 02: 6T 62 
oS vL Al | €L:80| € | 2M HS | 91:81 1 09| 2 TH St Og Fen) Ba Ss Zé | S¥:00 6 | W 6S: 61 Lé 
a 89 Aly |} ES2 ZL} +E | Som eS. | geri e or) = oT S = ri a ae ba U1 vy: €T 61 
I OV AI | Z€:80/| +€ | OPM IN | 00:80 £01090] 9 9 | OO:IL CO) au] S a £0 St | er: 00 SO 
= a — See oe) a 7 = a xs = =a zy — | penunuos s LG: 80 20 
= oF Al | ZT + 240 € | 824 GON | 00:20 82 | — ca = . - a ™ | peneyuao s 8¢ = SO 0€ 
= oT AI | ZV: 80 € | VIM GEN | SP:80 ¥ | — = <7 a = = —_ | penunuos SU | 9F:00 SZ 
9 , AE | #£°9% | +@-) OPH O2N | 80: 51-61} — = 7 00: UL ca) BE & SOT | ST: 2% 61 | — as 
Bu a 7 = = = 7 OD: | ql 9 00> GO GE} "a= os a a = a 
= GS AL} 22-02 | 6) S8a JEN?) OF SOL ot) — — = — +. a ed = s SE- LT ZT 
= OV AI | 10:00 £ | SEM 8IN | $0°00 8a) — aaa ed — a ak a = s vI: 91 10 
€I | OS at ace (€) | €€M ON | SO:91 12 | — = hae _ HE| 2 c9 | 00: S0 2 sm | 20:81 & 
cL | &@ Tae) coe OE & | SIM Sos | OL Tl O24 we - my . HL} 2 cL | 00: £2 02 s GE: er IZ 
ae GV AI | 81:2} @ | 2M 9IS | 0S: 12 se} — ar = > “iy Py = ur €0: ST 22 
ae ZS Mot SLO UE 1 E |) OSS SES: NEO MEL 08. ee . rr * = = a i ze 5} 0& - 0¢ 60 
cs 8E AI | 1&: éT S | Gls SIN. | ero is | tea tr co +? ou) S 6S | 00:91 I¢€ s 0€ : 20 20 
= 09 AI | €¥: dé @ | OH QIN | OC ize se) = “7 - — s = zy = ur 91: OL TE 
= 9S MI | 7e:ce| € | SIM HS | OE 2c ZI) = as = i. = wee saad — sw | 9:90 0¢ 
= 8S A) go2 28 | i OER Be | Oe SLES Lay fe _ = = — — rer SS s Jt = $0 Oi 
O'T ua AI | VS: 0¢ Go. | Peo SIN) OO 16.220) <= ai aie = aad 3 86 | 00:22 L2 Sw | 8&- 10 0€ 
= 0S AY.) Ug @ Geet. fa a Oaa|) = ahi = = —s a = a Sw | ah: €2 20 
S‘0 67 AI | S€:€0| +€ | S8M SZN | ¥E:€60 €2 | — _ 4 — HL; § = 1 OO) ar0cee S JO > 20. Sz 
na oe Al | #0:81|} +é@ | OLA 22N | SO: 8I 91 | — — = om; an = = a uw 1@: 0 61 
y y ty Bos | at mi? 9 Ot Ge NE Y Ue ee! ay Wy 
adAT| 442g | ‘dwy| uornrsog YRS * S90) ‘sqy | “ing wes «890 eddy | “ang cals epein qaes 
Vy Wad oIpey AIL] J Ia IWOLY | weis0Uu0y w401g 
sJUdAq IP[OS uoneziuoy deg arog oeuseUl0ey 


“unt 


“idy 
“IRI 
“Tey 
“Tey 


"qe 
8S61 


“AON 


20 
“das 
“das 
“das 
“das 
“das 


“ony 


“nl 
[nf 
“ne 
[nf 
‘unt 
“unt 
‘unt 
Ar 
aidy 
“Ie 
eek | 
‘uel 
LS61 
29q 
“AON 
‘das 
‘ony 
ACY 
ACI 
ady 
“IR 
“qo 
“qe 
9S6T 


ayeq 


‘anyL SAL ‘Aeg aynjosoy ‘ qy ‘molleg Julog ‘qd ‘euniry ‘yy "HORM FA hal ‘932110990 ‘OO x 


ror) 
low | 


OR N 
Ll Ileal 


BY: HPO Be 
loon lanowe | 
A 


K. SINNO 


2 


Se) 


iS ee 
Se oo en? Bh les | | | 


PERBEEE IEEE BPRPPREEERE PPRPPR IRERRRP Eee ERR ee | 


10 


T |:06M S2N |(ZT:02)02 | UM! S =~ (( 'sais [= — — — = — “AON 
+€ | SEM 9ZN | 40:20 ST | UM! ST — | 00:90 SI|/ am}| — | —-|,00:90 ST} #8 GG: 12 SI | “AON 
+€ | 2M OEN | 96:20 FT| — — ie — — — — — Sw | €0:€1 SI | “AON 
+€ | 2OM LON _| 22: E121 |} UM | OC! ap | vee SI 21 | uy] — =—'"| OO: Ot.2zt | 8 0€: OL €I | “AON 
(+8) |\(O@a O€N)|(ST: ¢0)IT | — = — — — _ — — s 8h: I ZI | “AON 

€ | 064 LIN | 0F:00 €0| — | — | — mie au| 9 09~| 0€:60 €0 | Sw /|(0E:7Z0)t0 | ‘das 

€ | SIM ON | S€:S90 #1} — | — ~ oe! au] s Ob | ST: EI #1 | SW | 60: FI OT | “Bny 
€ | €M ZN | OF: 12 zz} — | — | — — —_ — — — sw | g¢:61 62 | ‘unt 
€ | 904 O2N | 9E:11 Sz} — | — — — — — — _ Sw | 9p: 10 42 | ‘unt 
+€ | 9VA 82N | 62:80 ITO] — — — — a — _ _ sw | 0¢:2O 70 | ‘ung 
+Z | PIM ZIN | SS:80 92 |] — —|}— — — — —_ _ su | 61:01 82 | ARN 
+€ | SOM 62N | 22:S0 1] HL [s'P s9 | 02:90 €f| aN! 9D Sy A Mr 90 §I sm | 02:11 91 | AeW 
+€ | 800 O1S | *O:F1 90} HL | ST<| Tor] o¢:8t 90| an | — | or] — 90 | Ss 12:70 80 | AeW 
€ | 06M ZIN | ST:01 F0| HL! § vw | :01 1/4! S S GI: Il 70 | — ‘ae Ae 
€ | 2M TIN | €: 70 62) HL! PI vIT | 00:90 62 | auN!| §S os | 0€:80 62 | § €I: Zi 06 | “dy 
€ | wea SOS | 0€:10 82 | HL! ¢€ ve | 00:90 82|auN/ §S ST | 00:70 82 | Sw | Zé: 10 O¢ | ‘ady 
rd = = HL| & Ov | 00:80 SO | au |(9D)|} 2@E | oe f1 sO | w S€: 91 90 | “ady 
€ | 6OM SIN | S¥:80 10| HL! ¢€ 98 | Sv¥:60 10|/ au] S 29 | 00:01 10; Sm | €1:¢€2 Zo | “ady 
+2 | SIS TIN | OC:ST OF | HL Zz mI | 08:20 1€ | aM! Ss — Sh: @ Ik | 8 ev:iz2 I€ | “2eW 
€ | 1a ZIN | $0: 20 62] — |} — | — ne au| 9 LL | see rf an S$: 60 IE | “ze 
2 | 99M OZS | VEEL ST} — | — | — — — | — —~ — w Ge: 20 8I uef 
€ | 2a @N | 6:02 11} — |} — |] — - au| 9 LI | ces 91 SI | w 6S: 8T §I uef 
0961 

S | 2O0 GON | aiuto) — | + 7 = oo — — | — — — = ‘aq 
+I | 0664 PON | OL:61 10) — | — | — — —/}/—}]— — sm | 00:22 €0 | ‘das 

€ | MEM IIN | e2:01 81} — | — | =— — HL| 9 7 Ale Ol 61 | w 21:70 02 | “‘Bny 
+€ | 9M BON | ST? 12 91 | OD | SI<| v€ | o¢:22 9 | — |(S)| — | penunuos| s 86:91 LI "Inf 
+e | WA QIN | er: £0 FL | OD} St<! 18 | 00:40 $1 | — | — | — | penunuos! ¢ 20: 80 SI “Tne 
+E | OLA ZEN | OT:20 OT | O09] SI<| 06 | 00:4001|HL| 9 ail (laces 60 OT | sw | GZ: 91 II "(ne 
ae 06H LIN | 00:81 60} — | — | — — —}|}—}]— — sw | 61:60 IT unf 

974 22N | OL: SO €§1 
+€ | Za 80N | OL: 02 zt - Sum | €0:20 ST | AeW 
+€ | LPA €CN | 0S:02 OL | HL! SI<! O6f| O€:10 11 | ay! 9 091 | SI:ZO II | s 8Z:€2 IL | AeW 
£ | coy SIN | 00:2-21) — | — | = — —}/—|]— = w ev: Il #1 | “qeq 
6S6I1 
+2 | GOM €0S | zeéizi zr) — | = } = — —}—}] — — ul 22:02 SI | ‘20q 
€ | 9SM POS | Iv: FL FZ} — | — | — _ —-/}/—]}]— — Sur | €2: 6% | “Rg 
© | ZoM 20S) Ze reeiz] —} = | = — —/};—|]— _ su | 1€: 20 PO 
+@ | GM LIS | Th: 10 21 qal *# 89 | O€: FI 2 | HL|(9)] 22 | OO: 91 2 | sm | g0:%0 cz] ‘das 
+E [PLM 01S") 6:20 411 — | — a = — —-);—/— - w 0€ : 60 OT dag 
€ | VSM OZN | G0:00 92 | HL| st<} IZ | 00:1092!HL! S 0S | 00:20 92 | sw | €0:¢€0 22 | ‘BnYy 

S | SOM IGN. | JT eelrez AL) Se 08 | 00:21 2|auN| S 6S | 00:91 ce | sw | OF: 10 Fz | “Bny 

€ | 924 6IN | (EP00) 0¢ | da | ¢ OF 1.00 25Sr ie | ae |. 2 8I | Sk-ST I¢ | w 12:20 a | “sny 
+€ | €SM PIS | CE: 70 91 | HL | SI<X| 8% | 0€:90 91 | au | 9 69 | 0€:40 91 | sw | 72:90 ZI | ‘Bny 

€ | SVM FITS | £0: £0 62 | HL [20 2 | S0:'0 6¢/ au! Ss 42 | SI:%0 62 | w 62: ST Ig ny 
+€ | 2OM VEN | 6€:00 40] HL! SIX] 98 | 0€:20 20| aN] 9 €€I | SI:z0 20 | s 87:20 80 | ‘Inf 


The Solar Geophysical Events of November 1960 


By Tatsuzo OBAYASHI* 


Arctic Institute of North America, Washington D.C. 
(Read Feburuary 17, 1961; Received Aug. 1, 1961) 


Abstract 


A series of outstanding solar-geophysical events of November 1960 is des- 
cribed. Results obtained various measurments; solar phenomena, cosmic, rays, 
ionosphere, geomagnetism and aurorae are summarized and include a theoretical 
model of disturbances which appears to be consistent to date. It is shown that 
the events of November 1960 are unique in that they provide not only the in- 
formation of solar particles up to cosmic ray energy ranges, but also the knowledge 
of existing interplanetary magnetic flelds. Some important discoveries made in 


the events are also described briefly. 


A series of events occurred during November 1960 which were undoubtedly the 
most outstanding and complicated solar-terrestrial phenomena in recent years. Three 
solar cosmic ray events were observed which were associated with intense solar flares. 
The geomagnetic storms that followed were extremely severe. The upper atmosphere 
was very disturbed, and displayed unusually bright aurora and airglow. 

Major events which have been reported so far are listed in Tables 1 and 2. There 
were at least eight important solar flares which have been identified as the possible 
source of major terrestrial disturbances. All flares except on November 6, were ori- 
ginated from the same active region, the McMath Plage Region 5925, which passed 
the central meridian of the sun around November 12. Three of them, November 12, 
15, and 20, produced energetic solar cosmic rays, detected at sea level. Three major 
geomagnetic storms occurred that were accompanied by large Forbush type decreases 
of cosmic rays. Solar phenomena, geomagnetic activity, and records of cosmic rays 
and ionospheric absorptions (fmin) near the geomagnetic pole are illustrated in Fig. 1. 
Details of these phenomena are given in subsequent individual descriptions. 


1. The Event of November 12-13, 1960 


An intense flare started 13423m on November 12 near the central meridian of the 
sun, which had an active region a few days earlier. The flare developed from 13h25m, 
and reached a maximum at 1330m. An outburst of solar radio emission of spectral 
type II was observed between 13h27m to 13h31m at 200 Mc/s, and then a strong out- 
burst of type IV (continuum radiation) followed, reaching a maximum from 13h30m to 
13n50m, and fading out by 18h00m. The Ha intensity curve at McMath-Hulbert 


* Present address: Ionosphere Research Laboratory, Kyoto University, Kyoto. 
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Fig. 1. 


Table 1 Major Events, November 6-22, 1960 


* Time at 200 Mc/s. 
** Deep River 


; E 
Flare Radio: Outhorst SWF (SID) ie oe 
Date ] 
Time & : , * ‘ 
Imp position Type Time* Imp. Time 41% Time 
6 Il 1752-2030 | Maj., | 1827-1858 TY } 1708-1815 | = - 
O7E 13N 
10 Ill 1011-1430 - 1020-1116 I. | 1022-1152 - ~ 
28E 29N Maj., | 1116-1200 . 
11 I, | 0305-0428 Ill 0316-0330 Il, | 0311-0616 - (PCA) 
12E 29N II 0332 
IV 0340-0730 
12 Il, | 1323-1922 II 1327-1331 | IL, | 1325-1600 | 65 | 1340- 
4W 26N IV 1330-1800 | 120 1900-1030 
(13 th) 
14 II, | 0246-0520 ~ 0318-0335 Il 0300-0500 - (PCA) 
19W 27N IV 0335-0700 
15 Il, | 0207-0427 II 0221-0225 Il, | 0220-0630 85 0240-22 
33W 26N IV 0225-0700 
19 I 1543-1649 It 1559-1602 - - - (PCA) 
90W 28N (IV) 1636-1723 
20 Il 2017-2024 ; I 2023-2035 I_ | 2023-2145 5 2100-18.- 
(110W) 25N| IV | 2027-2046 | (21 st) 
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Table 2 Major Events, November 6-22, 1960 


han Geomagnetic Storm Associated Probable Delay 
Se Mag. Storm* Saas Pte Cosmic Ray Origin Solar Time 

Talis (4H)av. Time,,.i 4H )av, Forbush Decrease Flare hrs. 

y Til 

11 11-0033 12 1]---++s- 110 11-00::: M 6-1755 II 102.7 
12-13 12-1348 42 12-17-:- 450 bie = 10-1011 Il BG 
12-13 12-1846 36 12-19... 12-1930 Ss 11-0305 II, Boia 
13-14 13-1021 (200) 13-11--- (250) 13-1035 > 12-1323 IIl,. 21.0 
15-16 15-1304 22 15-17--- 240 15-1330 M 14-0246 II, 34.3 
15-16 | (15-2155) (10) 15-22... (15-2200) | (M) 15-0207 Il, 19.8 
21-22 21-0632 15 21-07--- 110 (21-06---) | (M) 19-1543 II 38.8 
21-22 21-2147 20 22-00:-- 70 (21-22---) | (M) 20-2017 II 25.5 


* Average of several low latitude stations, 4H in gamma. 


Ha Intensity 
McMath Hulburt 


100 


Radio Outburst 
Nera 2980 MC/S 


6000 Nera 545 MC/S 


1300 1400 1500 1600 


Nov. 12, 1960 
Fig. 2 Ha and radio outbursts records of the solar flare on 
November 12, 1960. (H.W. Dodson, and A. D. Fokker) 


Observatory and radio outbursts at Nera, Netherlands, are reproduced in Fig. 2. 


Tables 1 and 2: Information from Solar-geophysical Data, CRPL-F-B, NBS is included (Lincoln, 
1961). 
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As has been pointed out by Dodson, the Ha record of this flare resembles closely 
that of cosmic ray produced flares such as February 23, 1956 and July 16, 1959. In fact, 
the most outstanding feature of this flare was the production of energetic solar cosmic 
rays. At 13240m, about 10 minutes after the onset of the major radio outburst of aye 
IV, a distinct increase of the cosmic ray flux was observed at most neutron monitor 
stations in high latitudes. As shown in the record at Deep River, Canada, in Fig. 3, 
following a gradual rise, the flux intensity reached a maximum in about 2.5 hours and 
then decreased until a sudden rise occurred near 19h00m. The second increase had a 
maximum at 20h00m, exceeding 100 per cent above the normal level, and then a smooth 
recovery followed until about 10230m on November 13. 
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Fig. 3 Record of the standard neutron monitor at Deep River on 
November 12-13, 1960. (J. F. Steljes and H. Carmichael) 
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It has been noted that the solar flare and radio outburst intensities were sub- 
stantially reduced after 15h, and that there was no resurgence or new flare outbreak 
between 18h to 20h that could provide an obvious solar explanation for the second 
increase in cosmic rays. Therefore, the complex time variation of cosmic rays observed 
in this event must be understood in terms of the modulation effect upon the original 
flux of solar cosmic rays. 

In this respect, geomagnetic storms which occurred in this period are of particular 
importance, because they provide the information of any existing solar plasma cloud 
in space, and of the outer geomagnetic field, where solar cosmic ray particles might 
interact and likely be modulated. In fact, two sudden commencements (SC) of geo- 
magnetic storms were noticed in the earlier part of the event at 13h48m and 18h46m 
on Novemer 12. Magnetograms obtained at Hiraiso Radio Observatory are shown in 
Fig. 4. Unlike the first SC of the geomagnetic storm, the second one occurred right 
after the beginning of the main phase of the preceding storm. However, it was clearly 
seen in magnetograms obtained at most equatorial stations. Furthermore, a very sharp 
Forbush decrease of cosmic ray flux, which was observed in the meson monitor at 
MIT, followed soon after the second SC. Since almost all sudden Forbush decreases 
follow SC’s of geomagnetic storms within a few hours, it is rather convincing that the 
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Magnetograms Hiraiso Nov. 12 


Fig. 4 Magnetograms of H-component observed at Hiraiso Radio 
Observatory, geomagnetic latitude 26 N, on November. 12-13, 1960. 


SC at 18h46m is a really important one. It should be noted that the second large 
increase of cosmic ray flux coincides at nearly the time of onset of a geomagnetic 
storm. As will be discussed later, this fact is essential in order to explain the com- 
plicated structure of this solar cosmic ray event. 

As has been noted in Table 2, these SC’s probably originated from the earlier 
flares on November 10 and 11, respectively. Therefore, at the time of the flare of 
November 12, solar plasma clouds responsponsible for these storms were in the vicinities 
of the earth, but not yet enveloping it. It appears also reasonable to presume that a 
conspicuous sharp SC and a Forbush decrease at 10421m on November 13 were caused 
21 hours later by a plasma cloud ejected from the flare relevant to the solar cosmic 
ray event. 

In the polar ionosphere, soon after the solar flare of November 12, the polar cap 
absorption event started, which was noted by riometers as well as vertical ionosondes. 
There was also some indication that a weak polar cap absorption had already been 
in progress near the geomagnetic pole as early as 01h on November 12 (Gregory 1961). 
However, a more pronounced increase of absorption was certainly started around 14h 
00m. The riometer records at Cape Jones in Fig. 5, illustrate that the absorption at 
30 Mc/s went up steadily and reached beyond the dynamic range (12db) by 18h. The 
60 Mc/s riometers at Ottawa and Churchill, however, showed a sharp increase in 
absorption at approximately 1900m, coinciding with the time of the second cosmic 
ray increase. Most other stations also showed similar time variation, indicating again 
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Fig. 5 30 Mc/s riometer record at Cape Jones on November 12-17, 1960. 
(E. L. Vogan) 


a close association with the onset time of a geomagnetic storm (Vogan and Hartz, 1961). 

It is generally accepted that polar cap absorptions are caused by the precipitation 
of low-energy solar cosmic rays of 10-100 Mev. High altitude balloon observations 
made by Minnesota University group and rocket measurements at Churchill by NASA 
in this event, revealed the existence of an enormous amount of low-energy solar cosmic 
rays. The result from the Explorer VII (Van Allen, 1961) also confirmed this. The 
estimated peak proton flux above 30 Mev exceeded 10° cm~* sec™! and the exponent of 
integral spectrum in this energy range would be approximately of -2 to -3. 

Both polar cap and auroral absorption were extremely intense during the active 
phase of the geomagnetic storm of November 12-13, and a remarkable red aurora 
was seen as low as 40° geomagnetic latitude. The F-2 layer appeared to be completely 
wiped out during period of several hours, as indicated by ionosonde records and by 
lunar reflection experiments at Jodrell Bank and Sagamore Hill, Massachusetts. 


2. The Event of November 15, 1960 


A very similar, almost identical, intense flare was observed on November 15 at 02h 
07m. Radio outbursts started 02h21m which appeared at type II, and then almost 
immediately a gradual build up of a type IV continuum followed from 02h25m to 07h 
00m. The records of Ha intensity and of radio outbursts observed at Tokyo Astro- 
nomical Observatory and Hiraiso Radio Observatory are reproduced in Fig. 6. 

The cosmic ray increase associated with this flare (Fig. 7, at Deep River, Canada) 
started very sharply near 02h40m and rose, in 15 minutes, to a rather flat, jagged top. 
A smooth recovery began about 05200m at a rate similar to that of the event of 
November 12. Unlike the previous event, there was no SC of a geomagnetic storm 
until November 15 at 13h04m, and apparently the solar cosmic ray flux did not show 
any appreciable geomagnetic storm effect. 

On the other hand, the polar cap absorption observed by the riometer at Cape 
Jones (Fig. 5) showed that a very gradual rise after the flare continued until the time of 
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Nov. 15. 1960 
Fig. 6 Ha and radio outbursts records of the solar flare on November 
15, 1960. (S. Nagasawa and F. Yamashita) 
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Fig. 7 Record of the standard neutron monitor at Deep River on 
November 14-15, 1960. (J. F. Steljes and H. Carmichael) 


SC of the geomagnetic storm. At about 13h00m the absorption increased to very high 
values. Since there exists a marked sunrise effect in absorption, it is not certain 
whether the time variation of this polar cap absorption event well represents the 
change of solar cosmic ray flux in outer space, or the effect in the ionosphere itself. 
Though a geomagnetic storm started at 13h04m, the magnetic activity was fairly 
moderate until 21455m when a sudden outbreak of activity commenced (Fig. 8). This 
might have been the moment when the plasma cloud, emitted from the flare of 
November 15, arrived to the earth after a delay of about 20hours. This geomagnetic 
storm was rather small compared with the one of November 12-13, Therefore, auroral 
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20h46m. During this period, a remarkable ascending prominence was seen above the 
limb with a velocity of the order of a thousand km/sec. 

A distinct yet very small cosmic ray increase began at about 21h00m, which was 
a delay of about 30 minutes after the start of the radio outburst (Fig. 10). The rise 
time of approximately one hour, was comparatively slow, and after 23h00m, a gradual 
decay followed until around 18% on November 21. This is the first solar cosmic ray 
event whose source has been identified outside the visible solar disk. 
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Fig. 10 Records of neutron monitors at Sulphur Mt. and Deep River 
on November 20-21, 1960. (H. Carmichael, J. F. Steljes, D.C. Rose and 
B. G. Wilson) 


Theoretical Interpretations of the Events 


The most outstanding feature of the November events is the emission of energetic 
solar cosmic rays. Three of the events, November 12, 15, and 20, involved particles 
above the relativistic energies. There are some good reasons to believe that the flares 
on November 11, 14, and 19, produced the particles of sub-relativistic energies (low- 
energy solar cosmic rays). 

It has been known that most solar cosmic ray events are closely related to intense 
solar flares associated with type IV solar radio outbursts. In this respect, it is con- 
sistent that all type IV outbursts observed in this period were followed either by the 
cosmic ray increase of relativistic energies or by the polar cap absorptions, which is 
presumably produced by the particles of 10-100 Mev. Since type IV outbursts are believed 
to be the synchrotron radiation from relativistic electrons spiralling in magnetic fields, 
it is also likely that such agitated solar plasma bearing magnetic fields may be capable 
of accelerating solar protons from thermal to very high energies. However, it is not 
yet clear why only three of them produced relativistic solar cosmic rays. As has been 
stated by Dodson (1961), the flare of November 12 and 15 had very distinct features 
among other intense flares; they show the development, after flare maximum, of a very 
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complex loop-type prominence of absorption features and of large coverage of the 
nearby spots by a structureless extensive glow. The flare of November 15 was also 
seen by the naked eye in pearly white color, which might be of the synchrotron origin 
(Nagasawa et al, 1961). These facts are very encouraging for further studies of this 
problem. 

The time variation of cosmic rays on November 12 is unusual among other solar 
cosmic ray events in the past. The feature of double peaks is of particular interest 
because of its unique shape. As has already been pointed out, any obvious solar origin 
as the cause of this variation was ruled out since there was no resurgence or a second 
flare which might provide a new source for the second increase of cosmic rays. There- 
fore, the time variation of solar cosmic rays on November 12 must be understood by 
means of the modulation of flux in interplanetary space or in the outer geomagnetic 
field. 

Since a geomagnetic storm started soon after the flare, a depression of the geo- 
magnetic cut-off owing to the distortion of the field may cause an increase of incoming 
cosmic ray flux to the earth (Ionosphere and Radio Astronomy, Nera, P.T.T., 1961; 
Kodama and Kitamura, 1961). However, it is rather difficult to conceive the fact that 
a large second increase was also observed at Thule, near the geomagnetic pole where 
the effect is expected to be very scarce (Pomerantz, 1961). Hence, this mechanism may 
explain a part of the increase observed at latitudes outside the polar cap, but certainly 
not be responsible for the substantial part of the second cosmic ray increase. Thus, 
the modulation of cosmic ray flux must be taking place somewhere in outer space 
beyond the geomagnetic field. 

The argument follows that magnetic plasma clouds existing in the interplanetary 
space, play an impartant role in the modulation of solar cosmic rays. It has already 
been suggested by Gold (1959, 1961) that the solar plasma cloud ejected from the flare, 
may draw out any magnetic field existing in the vicinities of the sun, since the cloud 
itself has a high conductivity. As the plasma cloud advances into interplanetary space, 
it must form an expanding shell or a bulge of magnetic lines of force, such as shown 
in Fig. lla. Supposing that solar cosmic ray particles are injected into this magnetic 


a) b) 
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Fig 11 Schematic diagrams of the interplanetary magnetic fields : 
a) Magnetic cloud associated with a solar eruption, 
b) Remnants of magnetic fields, 
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bulge, then most particles are trapped inside by magnetic field lines, and advance in 
interplanetary space essentially with the same speed as the front of a cloud but not 
with the speed of the particle itself. When such a magnetic cloud evolves about the 
earth, which is indicated by the onset of a geomagnetic storm, a substantial rise of 
solar cosmic ray flux is expected. On the other hand, a magnetic cloud will also act 
equally to exclude the cosmic ray particles coming from outside the solar system, 
which may be identified as the so-called Forbush decrease of cosmic rays. 

The event of November 12 was exactly the same situation as described above. 
According to Steljes, Carmichael and McCracken (1961), the flare shot out solar cosmic 
ray particles when two magnetic clouds, which were presumably produced by earlier 
flares, were advancing in the vicinity of the earth but not yet enveloping it. Judging 
from the onset of two SC’s of geomagnetic storms, their distances from the earth at 
the time of particle injection were 10°km and 2X10’ km (1 A. U.=1.5x108km), respec- 
tively. Since the first SC was not followed by an appreciable Forbush decrease, it 
might be little magnetic field associated with it. However, the second one had a strong 
field so that it induced a large Forbush decrease. 

The majority of solar cosmic ray particles ejected from the sun were stopped by 
the inner (second) magnetic cloud, and were trapped. However, a part of the particles, 
presumably high energy ones, would leak out through the magnetic barrier and reach 
the earth before the arrival of the magnetic cloud. A slow rise and subsequent fall 
of observed cosmic ray flux can be explained by this mechanism. The maximum 
energy of particles produced the first part of the increase would be 5 Bev. The 
integral spectrum was undoubtedly very steep (approximately E~*). 

Immediately after the arrival of the second magnetic cloud, a very sharp rise of 
particle flux was attained. Although the rise was steep, the increase was simultaneously 
world-wide, and there was no impact zone effect because the particles had already 
well migrated from the entire magntic cloud. The large Forbush decrease, observed 
after 19h30m by meson monitors, certainly indicated that the earth entered within the 
magnetic cloud in which cosmic ray particles coming from the outside of the cloud 
were excluded. During the period of this second increase, the low energy solar cosmic 
rays became abundant. This was detected by rocket and balloon measurements. 

The time variation of polar cap absorption followed very similar pattern having 
two steps of the increase coinciding with the arrival of the magnetic clouds. This can 
be explained consistently by the trapping mechanism of solar particles, though the 
energy range relevant to this is of 10-100 Mev, considerably less than that detected 
by neutron monitors at sea level. 

The flare of November 12, which produced this cosmic ray event, also ejected its 
own magnetic plasma cloud. The cloud arrived at the earth 21 hours later, and induced 
an extremely sharp SC and a sudden Forbush decrease as well. Although the SC 
occurred during a very active phase of the previous geomagnetic storm, the sudden 
impulse was noticeably recognized all over the world. It was of considerable interest 
that the solar cosmic ray flux suddenly disappeared after the arrival of this plasma 
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cloud. It has deen found by Roederer et al (1961) that, by examining the energy spec- 
trum of incoming particles, the steep spectrunm characteristic for solar cosmic rays 
disappeared and a shallow spectrum remained which is common for galactic cosmic 
rays. 

This fact indicates further support of the existence of the third magnetic cloud. 
Solar cosmic ray particles remained in interplanetary space were swept out by this 
magnetic cloud. It is quite possible that in front of the sweeping magnetic cloud, 
solar particles may be piled up. This was found to be the case since the polar cap 
absorption showed a sharp increas just before the SC of the geomagnetic storm. 

On November 15, there was another increase of solar cosmic rays associated with 
an intense flare with a strong type IV radio outburst. Unlike the previous one, this 
had a much sharper increase and attained a single maximum within an hour. There 
was no geomagnetic storm until late on November, 15. However, at the time of the 
flare, the effect of Forbush decreases produced by previous magnetic clouds was not 
yet recovered. Thus, the earth was still within the huge magnetic cloud, and was 
connected directly to the sun by fairly well-ordered magnetic lines of force, as illustrated 
in Fig. 11b. Any particles ejected from the flare should be able to reach the earth 
quickly by spiralling down such lines of force. Also, their impact zones may be well- 
defined since the lines of force would direct and collimate the particles. In fact, 
McCracken (1961) has revealed that there was a discrete impact zone, and the particles 
were coming apparently from the direction of about 50° to the west of the sun-earth 
line. 

Nevertheless, there is a puzzling fact. An SC of geomagnetic storm occurred at 
13h04m on November 15. Judging from its delay-time, the plasma cloud responsible 
for this storm should have originated from the flare at 02h46m on November 14. If 
so, the plasma cloud was about 1/3 A.U. from the earth in the sun-earth line at the 
time when the solar cosmic ray particles were injected. Therefore, solar particles must 
have encountered this cloud before reaching the earth. The time variation of polar 
cap absorption did show some indication of this effect; a very slow rise of absorption 
until about the onset time of the geomagnetic storm. However, no appreciable effect 
was seen in the records of neutron monitors. It might be that the magnetic field in 
the plasme cloud was very weak, since the Forbush decrease associated with this 
geomagnetic storm was comparatively small. Therefore, the plasma cloud was rather 
transparent for cosmic ray particles, but still had enough magnetic fields to interact 
and trap low-energy solar cosmic ray particles. 

The November 20 event provides still another interesting fact. The ejection of 
solar cosmic ray particles has a fairly wide angular spread from the source, and an 
appreciable amount of particles can arrive beyond the visible solar disk. However, 
this statement may be valid only for the case of the flare in the western limb. Of 
course, statistically, it seems to be likely that solar cosmic ray events are larger for 
the source near the central meridian. 


It may also be worthwhile to note that the delay-time of geomagnetic storms, 
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which originated from the cosmic ray producing flares, was unusually short—21, 20 
and 26 hours, respectively, for the present three cases. Further, these geomagnetic 
storms occurred while the preceding storms were still in the very active phase. There 
is some statistical evidence that, when two plasma clouds are ejected from the sun 
one after another, the second cloud travels at a much faster speed than the first one 
(Obayashi, 1961). This might be related to some complicated situation in the inter- 
planetary magnetic fields, although a suitable explanation has not yet appeared. 

Several interesting discoveries related to the November events have also been 
reported. A strong enhancement of molecular nitrogen bands, especially the first 
negative band of N,*, 3914A had been detected in the polar cap region during the 
polar cap absorption events. Sandford (1961) has shown that the luminosity curve at 
3914A was very similar to that of the level of the polar cap absorption, essentially 
starting right after the solar flare. The variation was very steady, reaching a maximum 
approximately one-hundred times above normal at onset time of the geomagnetic 
storm, and then decaying to normal levels over the next two days. It appeared as 
an extensive glow rather than discrete auroral forms. The source of the emission 
exists at an altitude of 100km or below, and it seems likely that it is excited by 
incoming solar cosmic ray particles. 

Anomalous propagation of VLF radio waves has been confirmed during these solar 
cosmic ray events. It is believed that the radio waves propagate via an ionized layer 
formed by incident solar cosmic rays well below the ionosphere (Belrose and Ross, 
1961; Ortner, Egeland and Hultqvist, 1961). 

The material recovered from the Discoverer XVII, which was exposed to the solar 
cosmic ray event of November 12, has been under very careful study. An analysis by 
Yagoda and his group (1961) shows that there is a°considerable abundance of radio 
isotopes produced by the impact of energetic particles, which is presumably by solar 
cosmic ray particles. A large amount of tritium content was also found. Since the 
observed ratio of tritium to argon 37 is much larger than that expected from the 
experiment of particle bombardments, Fireman (1961) has suggested that production 
of the large tritium content is not from the bombardment by solar protons, but from 
tritium in the flare itself. 

Severe geomagnetic disturbances during the November events have offered the 
opportunity of investigating the geomagnetic storm effect on the air density of the 
upper atmosphere as revealed by the orbital motion of satellites. Jacchia (1961) and 
Grove (1961) found that appreciable increases of air density at the 200 km level occurred 
on November 13 and 15-16, which coincided with two maxima of the geomagnetic 
activity. 

Several other prominent features in the November events are strong absorptions 
in VHF ionospheric scatter propagation in the polar region, hydrogen emission in 
auroral spectra, and formation of remarkable auroral red arcs. Such phenomena seem 


to be associated with unusual severe disturbances during this period. 
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Abstract 


The entry of high energy charged particles is studied taking into account the 
effect of field distortion due to a geomagnetic storm. By using a modified version 
of Stérmer’s theory, the cut-off rigidities of incoming particles are computed for 
two models of the outer geomagnetic field. The result is compared with the 
observed effect on precipitation of solar and galactic cosmic rays during the initial 


and main phases of a geomagntic storm. 


Introduction 


Charged particles entering from outside the earth are always affected by the 
outer geomagnetic field. Recently it has become evident that the upper atmosphere 
is exposed to various solar and cosmic particles of a wide energy spectrum ranging 
from a few Kev. up to relativistic energies. The precipitation of low-energy solar 
cosmic rays has been investigated by high altitude balloon observations (Winckler 
and Bhavsar, 1960; Winckler, Bhavsar and Peterson, 1961), and by analyzing world- 
wide ionogram data (Obayashi and Hakura, 1960). The results have shown that the 
variations of solar cosmic ray flux, in time as wll as in space, are very sensitive to 
distortion by the outer geomagnetic field, particularly during the main phases of 
geomagnetic storms. It has also been revealed that galactic cosmic rays tend to increase 
at low latitudes during a severe geomagnetic storm, although the time variation is usually 
superposed on the Forbush type decrease (Kondo, Nagashima, Yoshida and Wada, 
1960). 

A tentative calculation of the particle precipitation into the distorted geomagnetic 
field was made (Obayashi, 1959; Obayashi and Hakura, 1961). It has been concluded 
that the shift of polar blackouts towards the equator as well as the storm-time cosmic 
ray increases during geomagnetic storms are explained consistently by the lowering 
of the geomagnetic cut-off for incoming particles. The present paper is an extension 
of the previous study and emphasizes that the theory can account for the effect of particle 


precipitation, not only in the main phase, but also at the initial phase of a geomagnetic 
storm. 


Theoretical Considerations 


Following the previous paper (Obayashi and Hakura, 1961), the geomagnetic field 
is postulated as having a potential combining a dipole field of moment M and a 


( 26 ) 
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uniform field 4H coinciding with the direction of the dipole axis. Using polar coordi- 
nates (r, 0, ¢), the field components H, and H, are expressd as 


H,=—(2M1— a1)cos a 
and H,=—("+4H)sind (2) 


During a geomagnetic storm, the impresse field 4H is positive at the initial phase 
and negative in the main phase. It is further assumed that the geomagnetic field is 
confined within a cavity of a certain size. The geomagnetic field beyond this cavity is 
regarded as a disordered field such that it has no systematic influence upon incoming 
charged particles. For the initial phase of a storm, the radius of the cavity R, is 


defined as the distance where Hr=0, viz., R=(S5) For the main phase, R, is 


defined as the distance where H,=0, viz. R= (40). Two models of a magnetic 


field are shown in Figure 1; (a) represents compression of the geomagnetic field, 
since there is no line of force across the surface r=R,, and (b) shows the effect of 
expansion of the geomagnetic field, there by reducing the field intensity near the 
geomagnetic equator. 


Fig. 1 Two models of distorted geomagnetic fields 
a) Model for the initial phase of a geomagnetic storm, 
b) Model for the main phase of a geomagnetic storm. 


Using these magnetic field models, the cut-off rigidities of incoming particles at 
various geomagnetic latitudes are computed. The cut-off rigidity of incoming particles 
P, at the geomagnetic latitude 9 is given by 


M Pte Bs a 
ae ES pe 3 
P, at (1 — cos‘@ Ca) 


where a is the radius of the earth, H,=0.31 gauss, and y is an angular momentum 
constant, being a complex function of the magnetic field intensity and the particle 
rigidity. 

Results of the computations are shown in Figure 2, where the relation of 4H and 
® for various cut-off energies of incoming protons is given. It is evident that for the 
enhanced 4H——the initial phase of a storm, the zone of particle precipitation shifts 
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Fig, 2 Geomagntic cut-off energy of protons of various latitudes for 4H, the 
Dst field, between +1007 and —400;. 


towards high latitudes; while for the depressed 4H——the main phase, it shifts towards 


low latitudes. The latitude shift is more pronounced for the lower energy particles 
than for the higher energy ones. For very low energy particles, 


1/3 


3 | 4H 


cos°O= : 
2} Ay 


This relation is valid only for the depressed field, and 4H may be regarded as the 
Dst-field, the storm-time variation of geomagnetic storms (Chapman and Bartels, 1940). 

Since solar cosmic ray particles have a very steep energy spectrum, a sharp latitude 
cut-off exists. Figure 2 is, therefore, useful for the study of such solar cosmic rays, 
and it is possible to predict a change of the precipitation zone from the observed, Dst- 
field of a geomagnetic storm. 

For the galactic cosmic rays, which have a shallow energy spectrum, the diagram 
shown in Figure 3 is appropriate. The change of cut-off rigidity 4P/P versus 4H is 
given for different latitudes. Normal cut-off rigidities (for the centered dipole field) 
are also shown. To determine the change in cosmic ray intensity from the estimated 
4P/P, one may utilize the observed latitude curves of the normal cosmic ray intensity, 
or else the known rigidity spectrum of the incoming particles. It is shown that for 
enhanced 4H——the initial phase of a storm, the cut off rigidity increases and con- 
sequently the cosmic ray intensity decreases; while for the depressed 4H——the main 
phase, the lowering of the cut-off rigidity results in an increase of cosmic ray intensity. 
At high latitudes, however, the effect of a geomagnetic storm would be very small, 


since there is no appreciable latitude variation of the cosmic ray intensity above 
geomagnetic latitude 50°. 


Discussion of the Results 


Particles precipitating upon the geomagnetic field behave differently according to 
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Fig. 3 Change of cut-off rigidity 4P/P for various 4H at different geomagnetic latitudes. 


their energy density. It is important to know whether incoming particles behave 
essentially like a single particle or like a conductive plasma. As has been shown by 
Ferraro (1952), below a certain critical particle density, the electromagnetic interaction 
among particles becomes negligible, and their motions may be treated as for a single 
particle. This criterion can be obtained approximately by equating the shielding depth 
of a plasma A, and the impact parameter of a particle in the geomagnetic field, i.e., 
p=, [me Ee (4) 
8xe?N Nes 
where J is the particle density and e is the charge in e.s.u. It has been shown else- 
where (Obayashi and Hakura, 1960) that the flux densities of galactic cosmic rays as 
well as solar cosmic rays above 10-100 Mev are sufficiently low permitting their treat- 
ment as single particles. On the other hand, the cloud of low energy particles ejected 
from the sun behaves like a conductive plasma, and therefore, produces magnetic 
distortions when it invades the geomagnetic field. Thus, the interesting situation 
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Fig. 4 An example of the geomagnetic storm effect 
on polar cap absorptions (ionospheric fmin) and on 
galactic cosmic rays. 


prevails that although a low energy solar plasma interacts with the geomagnetic field 
producing a geomagnetic storm, high energy particles penetrate through the distorted 
geomagnetic field into the upper atmosphere. 

A typical example of such a geomagnetic storm effect, for incoming high energy 
particles, is demonstrated in Figure 4. An intense solar flare, at 0830m on July 3, 1957, 
produced both a polar cap absorption, that is a solar proton event, and a geomagnetic 
storm. Enhanced ionization in the polar ionosphere started right after the flare, and 
the sudden commencement of a storm (SC) occurred 40 hours later at 00h42m, on July 
5. The records are reproduced of the H-component of the geomagnetic field at 
Honolulu, and of the ionospheric absorptions indicated by fmin observed at 14 high 
latitude stations in the northern hemisphere. 

The fmin data are divided into two degrees; blackouts, indicated by thick line, 
and The fmin above 3 Mc/s, indicated by thin lines. The polar cap blackout developed 
before the onset of a geomagnetic storm with the southern-most extent approximately 62° 
N. At the time of the SC, the absorption was reduced suddenly at latitudes below 75° N, 
and continued until the beginning of the main phase of the storm at 03230m on July 5. 
The southern-most extent of increased absorption was as low as 58° N during the 
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main phase of the storm. 


Since polar cap blackouts are caused by energetic solar protons precipitating into 
the polar ionosphere, the observed result is exactly what has been predicted in the 
present theory. The estimated boundary of precipitating particles of 50 Mev is shown 
by a dotted line Figure 4. Good agreement exists between the theory and the obser- 
vational data. Here it has been assumed that even before a geomagnetic storm, the 
outer geomagnetic field is distorted approximately 4H=—30y; by interaction between 
the geomagnetic field and the interplanetary gas, due to the continual orbital motion 
of the earth (Obayashi, 1959). 

A moderate Forbush type decrease of galactic cosmic rays was associated with 
this geomagnetic storm. The maximum decrease of cosmic ray intensity at high latitudes 
was approximately 4 per cent. At low latitudes, the decrease was not pronounced, 
but the intensity appeared to increase slightly during the main phase of the storm. 
This is also what was expected in the present theory. 

In order to see the effect-due to geomagnetic distortion, an intensity difference 
between Mt. Norikura and Resolute Bay is also shown in Figure 4. Since Resolute 
Bay is situated in the polar cap region where the change of geomagnetic cut-off rigidity 
barely affects the intensity, this difference represents the net effect due to distortion 
of the geomagnetic field. It is clear that the galactic cosmic ray intensity decreased 
slightly at the initial phase, and then increased appreciably during the main phase. 
Assuming the rigidity spectrum at low latitudes as P-'!*, the change in cosmic ray 
intensity is computed by 4//J=—1.5x4P/P,. For the initial phase, 4H~+50;, the 
expected decrease is 1.5 percent at Norikura (0=25.6° N). At the main phase, 4H=— 
1007, the increase is approximately 4 per cent, indicating a good agreement with the 
observed result. Therefore, it seems to be very probable that the storm-time variation 
of cosmic rays at different phases of a geomagnetic storm can also be explained by 
the effect of the induced Dst-field. 

Thus, independently observed facts concerning the effects of solar cosmic rays and 
galactic cosmic rays substantiate the present theory. However, some effects cannot 
be explained by this theory. Occasionally, polar cap absorption events show a definite 
expansion of the precipitation zone toward the equator at the onset time of geomagnetic 
storms. It is likely that in these events, energetic particles are trapped in the solar 
plasma cloud and consequently there is an enormous sudden increase of incoming solar 
protons upon arrival of the cloud. In such cases, the effect of geomagnetic distortion 
would be completely masked. Brown et al. (1961) observed a burst of X-rays along 
the auroral zone at the SC of a geomagnetic storm. They explained this pheno- 
menon as a dumping of electrons from the outer radiation belt. Winckler and his 
group (1960, 1961) claim that the effect of geomagnetic distortion for solar cosmic rays, 
viz., a lowering of the cut-off rigidity, is appreciable only during the earlier part of 
the main phase of a geomagnetic storm. Kellogg and Winckler (1961) has given a 
tentative explanation of this as the effect of inward shift of the geomagnetic ring current 


after the main phase of a storm. 
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Abstract 


Since the discovery of the cosmic-ray storm by Forbush in 1937, the various 
attempts have been made to explore its nature. Yet its very cause seems to be 
unclear in view of the lack of phenomena which are intimately correlated with 
the cosmic-ray storm. In this study the cause of the cosmic-ray storm has been 
sought after from observed evidences. 

In order to find the cause of cosmic-ray storms, the correlation of the spectral 
type of solar radio outbursts and cosmic-ray variations and magnetic activities 
has been studied for the data during the period from July 1957 to December 1960, 
The results obtained are as follows: 

(1) Cosmic-ray storms are closely associated with the radio outbursts of Type 
IV (continuum). The result indicates that the eruption followed by the radio out- 
burst of Type IV is the cause of cosmic-ray storm. 

(2) The size of cosmic-ray and magnetic storm was compared with the 
heriographic longitude of the eruptions. Large magnetic storms occur after 
eruptions near the central meridian, while the size of cosmic-ray storm is inde- 
pendent of the heliographic longitude of the eruption. 

(3) From these results, a possible model of magnetized corpuscular cloud 
is presented which seems to be ejected from the sun with eruption followed by 
Type IV outburst. This model of magnetic cloud is a large total bulk with a 
core and appears to be able to consistently explain not only the facts above 


described but also the other earth storms. 


Introduction 


A world-wide decrease of the cosmic-ray intensity accompanied with a magnetic 
storm was first discovered by Forbush“ in 1937, and this phenomenon was called the 
cosmic-ray storm. However, the correlation between the cosmic-ray and the magnetic 
storm is rather complicated; large cosmic-ray storms are often associated with very 
small magnetic storms, while magnetic storms occur without association with cosmic-ray 
storms. Therefore, it will be difficult to consider that the magnetic storm is a direct 
cause of the cosmic-ray storm. 

In 1953, an interesting statistical study was made by Sekido et al“, They classi- 
fied the magnetic storms into two different types, S-storms and M-storms, respectively ; 
the S-storms were closely related to sunspot groupes, while the M-storms were not, 


( 33 ) 
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They concluded that the former could be attributed to the magnetic corpuscular cloud 
ejected by solar flares and the latter caused by the non-magnetic corpuscular beams 
continuously emitted from the solar M-region. This implies that the cosmic-ray storm 
is caused by something related to sunspots, because magnetic fields should be responsible 
for modulating cosmic-rays. Actually, however, large flares are not always accompanied 
by the S-storms. Thus, it is needed to specify the flares by any other factor and to 
find out the individual correspondence, not statistical, between flares and cosmic-ray 
storms. 

In association with solar flares, on the other hand, there occur outstanding radio 
outbursts. Paying attention to this fact, we choose those flares which are responsible 
for the cosmic-ray storms, that is, the flares associated with Type IV radio outbursts, 
as described in the preceding paper“. In fact, a remarkable correlation has been 
found between cosmic-ray storms and the flares associated with Type IV (continuum) 
radio outbursts, but the size of the cosmic-ray storm is not always correlated with 
that of the magnetic storm. This seems to require another parameter which charac- 
terizes this correlation. 

The individual correlation between cosmic-ray storms and Type IV bursts makes 
it possible to locate the flare responsible for a cosmic-ray storm. The location of a 
flare may be a useful parameter which explains the relation between the sizes of the 
cosmic-ray storm and the magnetic storm. 

In the present study, we try (1) to confirm the correlation between cosmic-ray 
storms and Type IV solar radio outbursts, (2) to find any relation of the sizes of cosmic- 
ray storms and magnetic storms to the location of responsible flares, (3) to construct, 
in the light of the evidence presented, a model of the magnetized cloud supposed to 
be ejected from the sun. 

The data used in this paper are as follows: (1) Cosmic-ray data were obtained 
from the neutron monitor at Mt. Norikura. (2) Solar flare, outburst data were obtaind 
from Tokyo Astronomical Observatory, Hiraiso Radio Wave Observatory and Harvard 


spectral data. (3) Geomagnetic data are Kp-index and records at Kakioka Magnetic 
Observatory. 


1, Simultaneous Occurrence of Cosmic-Ray Storm and Type IV Radio Outburst 


In the previous paper“, we have pointed out that during the period of IGY, almost 
all cosmic-ray storms are associated with the Type IV outbursts. Further examination 
of available data shown in Table I. confirms this suggestion. Outstanding Type IV 
bursts were selected from the reports of references (6) (7) and (8), and cosmic-ray 
storms were defined as variations in which intensity decreases (47) relative to the 27- 
day running average value were 1 % or more. Magnetic storms were chosen so that 
they were preceded by Type IV outbursts. Picking up only such cases that two suc- 
cessive storms were well separated in time, we found that in 52 cases among 58 
Type IV radio outbursts were followed by the cosmic-ray storms during the period 
from July 1, 1957 through December 31, 1960. Likewise only of 9 cases among 61 
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Table 1 
Date Type IV outburst Flare oo span Reference 
Start Duration Imp. Position Imp. 1% H 
hm m 
1957 Jun. 28 | 0716 70 3 N09, E27 2+ 3.0 241 
Jul. 3 | 0849 25 3 N09, W42 2+ * 155 
16 | 1801 35 | $29, W30 2 0.5 90 
24 | 1802 73 3 524, W2a..* 8 On 52 
Aug. 3 | 1720 cont. 3 N25, E26 1 2.8 62 
28 | 0930 75 s S30, E35).° 3+ 4.8 143 
31 | 1303 73 3+ N25, W022 3— * 200 
Sep. 2 | 1310 60 3+ N11, W26 1+ * 289 
11 | 0331 91 3+ N11, W03 3 0.4 486 
12 | 1515 241 3 Nll, W19 2+ * * 
13 | 1419 97 3 * * 
26 | 1948 cont. 2 Na E27) 63 2.8 
Oct. 20 | 1646 90 3+ $25, W40 3+ 5.4 50 
Nov. 24 | 0903 50 3+ S15, E370 3+ 4.6 183 
24 | 1811 32 3 S12, £23 1 * * 
Dec. 14 | 1238 17 3 N18, E77 2+ SAS 
17 | 0735 31 3+ N22, E44 2 * 
, 22.4 1712 56 S N18, W19 1+ * 49 
1958 Jan. 7 | 1908 66 3 N30, E48 2 0.8 
20 | 1453 cont. 2 N28, E28 2 1.6 
Feb. 9 | 2109 95 ~ 3+ Sloyn, Wit & 2 pe] 617 
Mar. 1 | 2041 13 3-— S17, W55 1 172 83 
23 | 1002 100 S $14, E77.-—3+ 3.6 54 
Apr. 25 1.6 40 
May 8 | 2227 cont. 3 N15, E22 1 1.5 
27 | 2252 cont. 3 No data 1.6 142 
Jun. 4 | 2140 36 3 N14, W58 1+ jv | 144 
5 | 0838 20 3 N27, W65 2+ * * 
5 | 1704 24 3 $21, E70 2+ * * 
6 | 0434 36 3 N15, W772 * * 
29 | 2025 cont. 2 N25, E80 1 1.8 
Jul. 7 | 0027 113 3+ N25, E07. 3 yee! 472 
29 | 0304 57 3+ Sid, W433 10 
Aug. 16 | 0438 57 3+ $14, W50 3 1.0 198 
19 | 2316 38 3 N17, E22 1+ * 74 
22 | 1440 135 3+ N18, W0O9 —3 2.4 120 
26 | 0019 130 3+ N20, W54"” 3 er 135 
Sep. 6 P 
F 14 | 0835 20 rs) Sile we 2 i A, 85 
Oct.. 2:1 1701 cont. 3 N179. E60, «2 iP 
ZL | Zoec 18 3+ Ss, Wez 2 0.8 270 
24 | 1442 31 3 S04, W58 13 a oe 
Nov. 8 : 
Dec. 12 | 1415 157 3 $01, WO8 2+ 14 132 
1959 Jan. 6 | 1715 105 cont. 3 1.0 102 
10 | 1414 298 cont. 3 * 
22 | 1547 13. _cont.1,2,3 2.6 114 
Feb. 9 | 1312 45 - N07, E84 2 a 176 
9 | 1400 188 cont. 3,3+ 
TZ | 2ous z 42 Niy, too. o $0 99 
Mar. 25 | 1556 484 cont. 1 ou 274 
Apr. 7 | 1505 220 cont. 1 2.0 151 
14 | 1824 Goat Feont. 1523 * 
21 | 1830 1 cont. 3 0.8 160 
May 10 | 2116 88 N19, E46 3+ 4.0 167 
11 | 2022 20 N08, E39 2+ 4.7 109 
Jun. 9 | 1714 46 cont. 2 0.8 62 
11 | 0120 1.3 
Jul. 9 | 2023 110 Ni4, E67 2 0.7 188 
10 | 2009 35 N22, E70 3+ * 
14 | 0338 —100 Ni, EOS: .o 6.8 533, 
16 | 2120 250 N17, W30 3+ 9.3 330 
Aug. 14 | 0140 65 N12, E28 2+ 2.9 
18 | 1029 120 N11, W38 3 * 107 
22 | 2140 225 2 * 
23 | 1400 605 3 * 
24 | 1400 640 3 * 
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Table I (continued) 


3 * (8) 
Rae ip wen | a | we [Bo 

Sep. 1 | 1914 36 120 A 148 @) 
19 | 1912 57 cont. 1,2,3 bh 8) 
Nov. 30 | 1750 110 eS 4 He 8) 
Dec. : 1400 >77 2,3 aa: 8) 
1960 Jan. : iS) 

170 1-3 N23, E03 1.8 
2 1653 iB 1-2 S10, W37 2.5 8 

1 2 [ 
ae 09 1356 20 2 Nu, E41 1.2 (8) 
Mar. 28 | 2051 190 2-3 Nil, E39 1.6 (8) 

29 | 0000 47 2.6 (8) (7) 

30 | 1526 454 1-3+ | N12, E13 2.8 (8) 
Apr. 29 | 0208 56 1+ , W20 1.3 
May 6| 1414 118 13 S10, E08 3.1 s) 
17 | 1755 34 1-2 0.3 ®) 
4 N18, W90 7 8 

1 | 2012 6 1 ; 
tak 25 | 1215 165 2-3 N22, E06 2.1 (8) 
25 | 1717 126 1-3 N19, WOl * (g) 
27 | 0018 31 a8 , E35 * (8) 
27 | 2150 44 3 N23, W23 2.3 8 

ul. 15 : 
hes 11 | 1929 9 1-2 N23, E27 1.0 (8) 
Sep. 3 | 0038 16 2 N20, E87 0.8 (8) 
4 | 0006 22 2 1.6 (8) 
16 | 1717 114 13+ | $21, E66 0.3 8 
oT 12 | 1753 6 N11, W24 1.2 (8) 
Nov. 12 | 1345 255 2-3 N27, WO 4.9 (8) 
19 | 1636 17 N28, W90 1.0 (8) 
19 | 1708 15 1-2 N28, W90 * (8) 
20 | 2027 19 2 N25, W90 + (8) 
Dec. 5 | 1834 24 3 N30, E90 1.3 (8) 
14 2.0 (8) 
26 1.3 (8) 


cosmic-ray storms were not preceded by the Type IV radio outbursts. Since Type IV 
outbursts and cosmic-ray storms were independently selected, it can be concluded that 
practically all cosmic-ray storms have nearly one-to-one correspondence with Type IV 
radio outbursts. The result indicates that the eruption followed by the outburst of 
Type IV is intimately connected with the cosmic-ray storm. 

On the other hand, it was shown, in the previous paper, that during the IGY, 
Type II (slow drift) outbursts were found to be not associated with cosmic-ray storms. 
To examine this fact, 69 Type II outbursts with importance 3 or 3+ were selected 
from the data observed on the Fort Davis records during the period of 1959-1960. From 
the superposed epoch method about cosmic-ray intensity in which the occurrence of 
the radio bursts were taken as the zero day, it is shown that there is no significant 
association between tne cosmic-ray storms and Type II outbursts. On the contrary, 
Thompson and Maxwell have shown that a statistical association should exist between 
Type II radio bursts and cosmic-ray storms. Since, however, the Type IV bursts are 
often preceded by the Type II bursts,“ the result obtained by these authors seems 
to be due to the occurrence of a Type IV outburst superposed on Type II. 
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2. Relation Between the Meridian Distance of the Solar Flares and the Sizes 
of the Cosmic-Ray and Magnetic Storms 


We first examine the heliographic distribution of the flares with Type IV radio 
outbursts. Figure I shows that those flares with Type IV outbursts followed respectively 
by cosmic-ray storms and magnetic storms are uniformly distributed over heliographic 
longitudes. 


NORTH 40° 


SOUTH 40° 
Fig. 1. The Heliographic distribution of flares followed by Type IV outbursts. 


In order to study the relation between the meridian distance of the solar flares 
and the sizes of the cosmic-ray storms and magnetic storms, cosmic-ray storms are 
classified into three groups according to their size, that is, great storms with |4/|=4%, 
medium storms with 4%>|4/|=2 % and small storms with |4/|<2 %, where 4I/ is the 
difference of the daily mean intensity from the 27-day running average. Magnetic 
storms are likewise classified into great storms, medium storms and small storms 
according to |4H|=3007, 300r>|4H|=1507 and |4H|<1507, respectively, where 4H is 
the maximum range of horizontal intensity during the storm at Kakioka. In order to 
avoid the ambiguity in deriving the size of decrease, storms having occurred in the 
recovering periods of preceding ones are excluded. Respective numbers belonging to 
these classes are listed in Table 2. 

The numbers of these events in respective categories against the meridian distances 
of the positions of corresponding flare are illustrated in Fig. 1 (a) and (b). It is clear 
that the size of the magnetic storm is strongly dependent on the meridian being re- 
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Table 2 
Great Medium Small 
Cosmic-ray 8 LT 27 
Magnetic Storm 6 14 22 


NUMBER 


MAGNETIC STORM 


MERIDIAN DISTANCE 


Fig. 2. Dependence of occurrences of 
cosmic-ray and magnetic storms against 
the meridian distance of the position of 
solar flares accompanied by Type IV out- 
burst. The size of cosmic-ray storm is 
classified as: black (=4 %) hatched (=2 % 
and white (<2 %) which are the differences 
of intensity from 27-day running averaged 
values, Magnetic storms are classified as: 
black (23007), hatched (=1507) and white 
(<1507) which are the maximum ranges 
in the storms. 


sponsible for larger magnetic storms. On the 
contrary, the size of the cosmic-ray storm is 
nearly independent of the position of the cor- 
responding flare. One might be surprised by 
the fact that the size of cosmic-ray storms 
caused by the flare occurring near the central 
meridian is not too large. The reason for this 
may be explained as follows: 

It has been shown by Yoshida and Wada“ 
that the world-wide increases of the cosmic- 
ray intensity occur associated with sever geo- 
magnetic storms, though they are usually super- 
It is 
called the storm-time increase. Kondoh et al“ 


posed on the Forbush type decreases. 


have given an explanation of this phenomenon 
as due to the decrease of the cut-off rigidity 
for incoming cosmic-rays caused by the geo- 
magnetic storm. When a solar flare occurs 
near the central meridian, a large magnetic 
storm takes place on the earth and a large 
Dst field thus produced results in a considerable 
increase of the cosmic-ray intensity at middle 
and low latitudes. Therefore, the size of the 
cosmic-ray storm appears to be outwardly 
small. 


3. Model of the Solar Magnetic Cloud 
Producing the Earth Storm 


By way of combining the evidences des- 
cribed above and other results obtained by 
many authors, a model of the magnetic cloud 
may be suggested. 


Firstly, the terrestrial disturbances caused by a solar flare with Type IV burst 


are the followings: 


(1) Cosmic-ray increase: solar particles of rigidity >1 BeV. 
(2) Sub Cosmic-ray increase: solar protons (~100 MeV) inpinging on the upper 


atmosphere at high latitudes.” 
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(3) lonospheric disturbance: polar cap blackout (~10 MeV), Type HI absorption“ 
(~10 MeV) and auroral zone blackout“) (~1 MeV). When a flare occures in the 
region of a high magnetic field of a sunspot group, an ionized gas cloud will be ejected 
outwards from the flare, in which a magnetic field will be frozen, so that high energy 
particles will be trapped. Type IV outburst is generally believed to be due to the 
synchrotron radiation from relativistic electrons trapped in this cloud at the corona 
region. Boischot“® has already proposed that during this stage, extremely high energy 
protons could be generated and then might escape along the magnetic line of force 
to give rise to a terrestrial cosmic-ray increase. This suggestion has been indeed 
justified even for lower energy protons by the above evidences (1), (2) and (3). 

Since it is accepted that cosmic-ray storms are caused by a magnetized cloud, the 
fact (3) is supposed to be a direct evidence that not only high energy particles but 
also magnetic clouds are ejected to the interplanetary space from the sun with Type 
IV burst. By Obayashi and Hakura the relationship among the various disturbances 
at the time of an earth storm has been classified, as shown in Table 3. 


Teble 3 (after T. Obayashi et al) 


Type of Combination of ) Magnetic Storm | Cosmic-Ray Storm Combined 
Terrestrial Event Type Type Type 
Type of Radio Outburst Ul IV IV 
Polar Cap Blackout Absence Presence Presence 
Auroral Zone Blackout Presence Absence Presence 
Magnetic Storm Typical Absence Presence 
Cosmic-Ray Storm Absence oe chelate 
Recent Event Nov. 6, 1957 Oct 2201950 Feb. 11, 1953 
Flare Importance 3+ 3+ 2 
Flare Longitude | E 27° W 40° Ww 14° 


According to this table, the authors have postulated that there are at least two 
types of solar corpuscular clouds which are responsible for the earth storm. The 
cosmic-ray storm type, characterized by the existence of high energy particles of 
10~100 Mev, is associated with the Type IV outburst, and the magnetic storm type, 
consisted of nearly thermal particles, is not, accompanied by the Type IV outburst. 
These typical two extreme types are, however, rather infrequent but the combined type 
actually takes place most frequently. 

Although these types of the magnetic cloud have been classified according to the 
conditions at the time of occurrence on the sun, the frequency of the high energy 
type seems to be too rare to explain the observed facts. Neverthless, it may be possible 
to explain this by taking account of the structure of the magnetic fleld that has a 
core in a large cloud. According to the facts described in the preceding section, if 
the core invade in the earth field, a rather large geomagnetic storm may be produced 
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and during its main phase both an auroral zone blackout and a storm-time increase 
in the cosmic-ray intensity may follow. Then the size of the cosmic-ray storm would 
not be so large. On the other hand, if the earth will break into the surroundings of 
the core, no appreciable magnetic storm may occur and consequently no storm-time 
increase is expected so that the size of cosmic-ray storm would be rather large. The 
occurrence of the cosmic-ray and sub cosmic-ray increases are independent of whether 
or not the earth penetrate into the core. 

Further details of the structure of the magnetic cloud may be inferred. The core 
has a considerably regular field which may be ejected from the sunspot field or their 
neighbourhood and in this field low and high energy corpuscules may be trapped. As 
far as the low energy corpuscules are concerned the core can be regarded as a dense 
plasma. The surroundings of the core has a irregular field, diffused from the central 
regular field, and high energy particles trapped in the core may gradually diffuse into 
the irregular field, as the cloud is travelling through the interplanetary space. There- 
fore, the particle density in the outer cloud may be very low. The geomagnetic storm 
would be produced by rather low energy corpuscules which concentrate in the central 
core and rather high energy corpuscules expanded in the outer cloud would give rise 
to the cosmic-ray as well as sub cosmic-ray increases. It seems to us, although not 
definite, what is essential to the cause of the geomagnetic storm is not the field intensity 
but the particle density, because as indicated by Sekido et al, even non-magnetic M- 
region would be able to produce magnetic storms. Hydromagnetic shock front preceding 
the magnetic cloud will be effective to the sudden commencement of the geomagnetic 
storm. 

It should be noticed that since all cosmic-ray storm are related to flares over the 
visible hemisphere, the dimension of a magnetic cloud may be of a considerably large 
scale ; its lateral dimension may be considered to be as large as one astronomical unit 
at the earth’s orbit. 


Conclusions 


In conclusion, the results described in this paper are as follows: 

(1) Almost all cosmic-ray storms are closely associated with the radio outbursts 
of Type IV. It is indicated from this that the eruption followed by the radio outburst 
of Type IV is the cause of cosmic-ray storm. Evidence has been obtained that the 
occurrence of a Type IV outburst is a useful event for the prediction of the geomagnetic 
storm and the cosmic-ray storm commencing 1-3 days later. 

(2) According to (1) the location of the eruption producing the cosmic-ray storm 
has been possible. Therefore, the size of cosmic-ray and magnetic storms were com- 
pared with the meridian distance of the eruptions. Large magnetic storms occur after 
eruptions near the central meridian, while the size of cosmic-ray storms is independent 
of the meridian distance of the eruption. 

(3) From the relationship between the location of the eruption producing the 
cosmic-ray storm and the size of cosmic-ray and magnetic storm, a possible model of 
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magnetic cloud has been presented. This model is supposed to explain a series of 
phonomena in the earth storm. 
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Abstract 


A theory of electromagnetic induction by a convectional fluid motion within 
the earth’s core is attempted in the hope of accounting for the localized intense 
secular variation in the geomagnetic field as has been found in the Antarctic 
area. Ifa toroidal magnetic field of 300 gauss at maximum is supposed to exist 
in the core, a steady convectional motion described by a spherical harmonic of 
order 5 and degree 5 can give rise to a poloidal magnetic field (degree 6, order 
5) of which the Gaussian coefficient of the magnetic potential amounts to 0.0457 
at the earth’s surface provided the radial velocity of the motion is taken as 
0.01 cm/sec. 

A study of the growth of the field tells us that a secular change of the order 
of 1007/yr is expected as long as a velocity of 0.1cm/sec is assumed. This 
order of velocity would not be impossible for such a localized motion as con- 


sidered here. 


1. Introduction 


Since Elsasser“? emphasized the importance of electromagnetic induction within 
moving material in the earth’s core, the possibility of accounting for the non-dipole 
part of the geomagnetic field and its secular variation through the induction process 
has drawn attention of a number of authors. For example, Bullard® has considered 
induction in eddies idealized as rotating spheres in order to account for the secular 
variation observed in South Africa and its neighbourhood. In his result has been 
suggested an important conclusion that a toroidal field much stronger than the dipole 
one is likely to be in existence in the core. 

Now since IGY, a fairly good network of observations of geomagnetic secular vari- 
ation has been extended over the whole surface of the earth, and various features of 
geomagnetic secular variation, especially in the southern hemisphere where only very 
little had been known until that time, have steadily been clarified. As pointed out 
and confirmed by one of the present writers@»“.@), a remarkable fact obtained in the 


recent studies on the secular variation is that an anomalously large secular variation 


( 42 ) 
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amounting to about 200 7/year in rate of change is taking place in the East Antarctic 
area. 

It has been concluded that this large secular variation has continued at least 
during the past 20 years. The isoporic chart of Z for the period of 1955-60 in the 
southern hemisphere is reproduced in Fig. 1, where the magnitudes of X and Y are 
also illustrated by vector arrows. 


Fig. 1. Isoporic chart of Z for 1955-60 in unit of r/year 
in the southern hemisphere. Full circles show observing 
points and arrows represent H=(X, Y) at those points. 


180° 


Fig. 2. Equivalent current system of the geomagnetic 
secular variation at the depth of 2900km, in 
unit of 10* Ampere/year. 
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The isoporic chart being subject to a spherical harmonic analysis of 6 in its highest 
degree, the equivalent current system of the variation assumed at the depth of 2900 km 
(namely, on the surface of the earth’s core) is obtained as illustrated in Fig. 2. As 
shown in this figure and as confirmed by numerical values of the spherical harmonic 
coefficients, the localized distribution of the intense secular variation in the Antarctic 
area can roughly be represented by harmonics of 4, 5 and 6 in their degree. 

Then, the question of how such a localized intense secular variation can theore- 
tically be explained may rise again. As a cause of secular variation, the effect of 
the westward drift of non-dipole components of the geomagnetic field has widely been 
accepted“), Provided that the origin and a sufficiently long life of non-dipole magnetic 
fields produced in the earth’s core are reasonably explained, a possibility of the west- 
ward drift of these fields can be theoretically understood™. According to the result 
of analysis, a tendency of westward drift of 0.18 degree/year in average velocity (ranging 
from 0 to 0.6 degree/year) can be detected for harmonics of »=2~4, but not for 
harmonics higher than 4 in their degree. Further, comparing the localized secular 
variation directly to the non-dipole field® in the Antarctic area, it is found that the 
hypothesis of westward drift of the non-dipole field can not account for the cause of 
the secular variation. 

Then, the most probable theoretical explanation of this phenomenon may be to 
assume upwelling of a part of strong toroidal magnetic field caused by a particular 
fluid motion within the earth’s core, as already suggested“. 

Bullard’s work® for interpreting the South African secular variation was along 
this direction of research. The induction process by rigid rotators embedded in station- 
ary conductors, such as assumed in the Bullard’s first paper, has been later examined 
extensively by Herzenberg and Lowes. In contrast to rigid rotator models, more 
realistic models based on interaction between magnetic field and fluid motion have 
been also studied. Coulomb“ has discussed the advection of the poloidal field on the 
basis of a simple model for upwelling motion in the core. The idea that a toroidal 
field which is presumably strong at some depth in the core can be brought up by 
convection has been discussed by Parker“. Although such a distorsion of a toroidal 
field by convectional motion seems to have an important bearing on the dynamo 
process which maintains the main geomagnetic field, no exact theory with detailed 
estimate has been so far available except the work made by Allan and Bullard@”, of 
which only a short note has been published. It is the intention of this paper to see 
what the extent of the poloidal field will be outside the core when a fluid motion of 
convection type and a toroidal field are specified in the core. And it is hoped to find 
some clue for explaining the intense variation in the Antarctic area. 


2. Theory 
The induction equation for magnetic field h is given by 


dh/Ot—(4ne)"'p2h=curl(v A Hp) beta) 
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where o, D and Hh denote the electrical conductivity, velocity of the fiuid motion and 
stationary magnetic field. It is not intended here to take the equation of motion into 
account, so that the present treatment is not really hydromagnetic. Some discussion 
of the hydromagnetic effect will be made later. 

The fluid motion is assumed as a poloidal one. The 7, 0 and ¢ components of 
the velocity are given as 


ns Rae tye 
dé% m n- EOYs 
| +(m+1)&% |r “aa 


den oy™ 
— 1)eé" n-1 OL n _ 
—|r& +(m+1)&? |r sind6¢ 


v ey) 


where Y% stands for either P?(cos@)cosm¢é or P%(cos0)sinm¢é and the origin of the 
polar coordinate is taken at the centre of the core. 

‘The theory of the origin of the earth’s magnetic field suggests that there is a 
very strong toroidal field, presumably as strong as several hundred gauss, in the core. 
Let us then assume that the stationary field is given as 


0 
z—0 (3) 
Tyr) 


The induction process due to D and Hi that is shown by the righthand-side of 
equation (1), then gives rise to magnetic fields of various types. 
They are 


nw 


Poloidal : ens Lard 

Toroidal : —29 a i hugl es 
where m and m denote respectively the degree and order of spherical surface har- 
monics involved. The fields of toroidal type do not appear outside the core provided 
the conductivity in the mantle is assumed to be very poor, so that our concern is 


confined to the poloidal ones which are readily written as 


Sey field : 
Rose Pha Y ens 
24 ast m | n— mee 1 
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where 
Yu=ayn/d¢ (6) 


Introducing (4) and (5) to (1) of which the righthand-side has been calculated from 


(2) and (3), we obtain differential equations for radial parts sp-1 And S24 ads 


d?sm-1 , 2M ASn-1__pogm  — n— es) het 
dr +r r ar . pe ; 
d?sms1 4 nt2) dss _ Rsta1=Snril”) ited 
dr° r dr 
where 
k?=4rop : : 


and p is algebraic operator 0/0t. fr-, and fry: are derived from the induction effect 


and given as 


fra(=— gn ai" Anaé™(nTAnr? (10) 
= = Se n n—m+1 m 14 
Sarl?) eee er Anoén(r) T(r) (11) 


The solution of differential equation (7), that is finite at r=0, is obtained as 
smei(M= ry [ Ans nth) t Inari errr K yaar 
niahr)\ "Dy lk fn(r)dr (12) 


in which A,-, is a free constant which is to be determined by the boundary conditions. 
J and K are modified Bessel functions. 

If we ignore the electrical conductivity in the mantle, the potential of the magnetic 
field corresponding to (4) is given as 


we-.=a(r/a)-4 YS (13 ) 


outside the core, where a denotes the core’s radius. i”, is then the coefficient for 


the magnetic potential originating from the interior of the core. The magnetic field 
is then obtained as 


eet 
ie tril m1 


ht_,= (z ae il oe (14) 
. i Ne fa oY ees 


At r=a, the continuity condition for the normal and tangential components of 


the magnetic field (the magnetic permeability is assumed as unity in electromagnetic 
unit) then leads to 
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—(n—1)ns%_(a)a"?=nir_, (15) 


dst m™ n— “me 
(78 ¢nsts) a Sat ( 16 , 
from which, eliminating i”_,, we obtain 
dsm; ] 
r——* + (2n—1) sm == () 
[ 7 +( )Smt | (17) 


This is the condition by which constant A,_; is to be determined. Before we 
proceed to determine A,_;, however, we have to perform the integrations involved in 
(12). Since no detailed distributions of €%(r) and 7;(v) have been so far known, it is 
necessary to make some assumption concerning these distributions. According to 
Bullard and Gellman“, the toroidal field in question is very small at the core-mantle 
boundary getting larger with the increase in depth. In one of their examples, the 
maximum of the field intensity occurs somewhere at the depth of 0.4a. Because of 
the solid inner core, however, it would seem more likely to suppose that the maximum 
occurs at a shallower depth. 

For the purpose of making the integrals in (12) tractable, it is assumed here that 
the toroidal field is zero except in a spherical shell between r=c and r=c+dc, dc 
being small compared to c. c=0.8a@ and 4c=0.1a@ as chosen in the later numerical work 
would be one of the likely combinations for actual values. In that case, (12) reduces to 


Sa-(n=r "t2{ Ag sl 1 (kr) +0" 2 Ac fas {In—1/R1) Kn—12(Re) 
— Ky RY) In-1/2(ke)} | (18 ) 


where f,_; is the mean value of f,-; throughout the shell being approximated by 


at ——3%t1 nt+m , E™¢)T»(c)c? 19 
F nt ‘sea, a WG Paiva (19) 


Going back to condition (17), we may now determine A,_,. We obtain 


A,i= —e Act, 4 [Ky-1 2c) +] n—1/(ke)Kn—s/(Ra@)/In—s/2ka)] (20) 


m 


Putting (20) into (18), s%_,(7) is finally solved. In the following, only sy_, at r=a 


is written down as 


Cc n—1/2 = 1 i fee. / (Re) 
m ny pte BOF 4 ad a ay 
Sn—i(@) ( =) 5 ch ; ka Tn-3/hka) 


With s”_,(a) thus determined, (15) enables us to obtain 77, by which we can estimate 
the S™_, magnetic field emerging through the core’s surface. 
As for the S”,, field, a theory similar to the previous one leads to 


c\ 73/2 au LiL as (kO) 
m (qy=—(& Ae Fy: oo 
Sn+i(@) (=) agi *Y ka Tn+1/2ka) 


where 


zg HM n—Mt! 47,8" ¢)T,(c) (23) 
Fast ee St an, : 


48 T. NAGATA and T. RIKITAKE 


3. Steady state 


Making p>0 or k>0, we have 


Pesnlle) Ty (ae) n—1/2 


Tn-3/2(ka)~ ie ) i " : 


Putting these into (21), we obtain 


Sn—1(a) = — (£)" ade fat for to (24) 


In the same way, we get 


c \ 2nt4 He Fi 
m pA pele adc! = for too (25 ) 
Snixi(@) =) 2n+3 
It is easily seen that solutions (24) and (25) can be directly obtained from differential 
equations for the steady state. 
Now we are in a position to estimate what velocity produces what magnetic field. 


Let us rewrite (23) as 


3 n—m+1 


7 Araés ath ig CFS ( 26) 
(n+1)(m+2) 2n4+1 


las 


where v,’ multiplied by Y™ is the r component of the velocity at r=c as given by 
ve =| nin +1) Er" | eam (27) 
r=¢ 


In numerical work, it is more convenient to make use of spherical functions defined 
by Schmidt instead of those functions of Neumann’s definition which have been adopted 
in the theory in the last section because, being partly normalized, Schmidt’s functions 
are always less than unity. €:,, in the above expression is introduced because of the 
partial normalization and is given by 


={tm 
Enym ie for m0 (28 ) 
Meanwhile, we take 

Ay! =[T2r*],<c (29) 


so that, assuming that the maximum value (at 6=45°) of the stationary toroidal field 
as 300 gauss, Hy’ is estimated as 200 gauss. 


On the other hand, the coefficient for the Stui field (Schmidt’s function is again 
used) is given by 


ee n—m--1 C\RE 
ln4i= 4 /H,! (£) n+lym 
ara (2n-+1)(2n-+3) ma v,' Hp! dc - —berlon (30) 


Enem 
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Suppose that there is a convectional motion of which the velocity is described by 
n=5 and m=5. (30) enables us to estimate 75. Let us take 


o=4X10- emu. 


4c=0.1la 
c=0.8a (31) 
a=3470 km 
In that case we see 
tet. (32 ) 
for 
v,’=0.01 cm/sec (33 ) 


The Gaussian coefficient of the potential at the earth’s surface is then given by 
(2)" ?§ (a): the earth’s radius) which amounts to 0.0457’. Judging from the observed 
magnitude of non-dipole fields, it is therefore possible to suppose that the velocity of 
the convection motion would be the order of 0.01 cm/sec which has been also suggested 
by Bullard and Gellman®® in their theory of geomagnetic dynamo. The induction 
process considered may well be responsible for the non-dipole fields. 


4. Non-steady state 


Going back to (21) or (22), which may be regarded as operational equations, we 
can examine the time-dependent behaviour of the S%, and S%,,; fields when some 
sudden change in the velocity takes place. If fn», changes instantaneously at t=0 
and remains constant afterwards, we obtain 


~ eee ae Fo 1 Insn(v/4nape) & 
a—1\4, L)=— | — Ac f na——\ a tN OE f f>0 (34 
Sn—1(4, 1) ey Pa =i bere Tn-213(/ 4napa) P P ST ge eee 


where L is Bromwich’s path of integration on the complex p-plane. 
It is seen that the integrand has a branch point at p=0 and simple poles at 


p=—a,?/(4nca’*) Ce oe ree (35 ) 
while a,'s are the roots of Jp-3/2(z)=0, so that 
Jn—a/2lAs)=9 0<a,<a,<a3< were ( 36 ) 


We take a path of integration as shown in Fig. 3 in which Bromwhich’s path of 
integration is shown by the straight line along the imaginary axis, so that all the 
singularities lie on the lefthand-side of the path. The branch point and poles are also 
shown by O, Pi, Po, Ps,.+0+0 . If we further consider the integration along the closed 


circuit like the one in Fig. 3, we get 


§2=Sot) p.t$ pt) pt sesees (37) 


where the integrals on the righthand-side are the ones integrated counterclockwise 
along the small circles centred at OP jes, PB tase . It is readily proved that the 
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Fig. 3. The path of integration. 


integrals along the large semi-circle vanishes when we make the radius infinite. Mean- 
while, the integrals along the two straight lines (parallel to the real axis) taken on a 
Riemann-plane are also proved to cancel one another. 

Putting (37), the integrals on its righthand-side being evaluated with the aid of 
the theory of residues, into (34), we obtain 


Sn—1(@, t)=Fy_ (t) st (a, ©) (38 ) 
where 
Fys()=1—4(2n—1)(£) 1 June) aie (32) 
a — In—1/2(@s) —Jn—s/2(@s) as 
In a similar fashion, s”,, (a, t) is obtained as 
Sn4i (a, t)= Fri (t) sti (a, co) (40 ) 
where 
Cc __Bs?t 
Pya)=1-4(2nt3(£)"" § Juvan( G6) @ te (41) 
a 1 Jn43/2(Bs)—Jn—1 PAR) ~ fe 
and 
Jn+1/(Bs)=0 0<Bi<Be<Bg<....s ( 42 ) 


For small values of ¢, (21) and (22) provide approximate equations like 
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e \ 1/2 = —k(a—c) 
Sri(a, p)=— (£) ee a 


Sil, == (Sede Fuys fre 
so that we obtain 
Fy1()=(2n—1)(£) me Berit —(1-£)(1-ers nba (a=c) 
—n—5/9 vie —c)? : ( 44 ) 
Fuss ()=(2n+3)(£) ef Z ghana da (1-£)(1-erf gaitir=e)) 


which are sometimes convenient for numerical work. 
Taking n=5, for example, the growth of the S$ field is calculated. F(t) as obtained 
from (41) and (44) is illustrated in Fig. 4. We see that the growth of the field is so 
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Fig. 4. Growth of the S,° field. 


slow in the beginning that no appreciable field is observed until 10? years past from 
the beginning. After that, the rate of increase gets larger, the 70 per cent of the 
steady field is attained after 4x10* years. After 10‘ years, the steady state is nearly 
reached. 


5. Discussion and conclusion 


A simple estimate of the time required for transportation of a toroidal field by a 
motion of fluid, the conductivity of which is so high that the magnetic lines of force 
are frozen in it, indicates that the time is given by d/v“” where d is the distance 
and v is the velocity. In order to bring up the toroidal field from a depth 0.8a to the 
core’s surface, it takes only 220 years for a velocity 0.01cm/sec. From what we have 
been dealing with in the last section, however, we see that it takes about 1400 years 
for the S field to reach 10 per cent of its steady value. It seems misleading to 
define the time for the growth of poloidal fields simply by d/v in this case. If trans- 
portation is made only by means of diffusion, it would take a time of order Arad?) to 
diffuse from a depth d through material of conductivity o. Taking o=3x10 e.m.u. 


as before, the time is about 7700 years. The induction process, therefore, seems to 


give a time for transportation shorter than the one in the case of simple diffusion 
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process. 
The curve in Fig. 4 indicates that the highest rate of increase for the radial 


component of the S$ field is of the order of 107/yr at the earth’s surface provided 
a S* motion having a velocity of the order of 0.01cm/sec is assumed in the core. 
If we assume 0.1 cm/sec for the radial velocity, secular change of the order of 100y/yr 
can be explained. According to the theory of geomagnetic dynamo“, maximum 
radial and azimuthal velocities have been estimated as 0.014cm/sec and 0.04 cm/sec 
respectively. For a more localized motion, it would not be impossible to suppose a 
velocity of the order of 0.1cm/sec, so that the induction model considered here may 
account for the extremely intense foci of secular variation amounting to 200 7/yr at 
maximum in the Antarctic region. 

The theory of induction has been in the above developed by specifying the fluid 
motion. Physically speaking, however, the driving force rather than fluid motion 
should have been specified because the pondermotive force produced by the growing 
magnetic field might become appreciable resulting in a considerable modification of 
the original motion. Such a dynamical consideration is a matter of great difficulty, 
no rigorous account based on exact magneto-hydrodynamics has been so far published. 
Even in this section, a very crude discussion of this point can be made in the following. 

Among the various terms of the equation of motion for the core’s fluid, it is known 
that the Coriolis force and the magneto-mechanical forces are the most important“ 
(8), For a nearly steady flow, the magneto-mechanical force is then approximately 
balanced by the Coriolis force, the magnitude of which is given as 2wpv (w: angular 
velocity of the earth’s rotation, : density of the core’s fluid). Taking v=0.1~0.01cm/sec, 
the Coriolis force and the magneto-mechanical force working on a unit volume of 
fluid are of order 10~-* ~ 10-5 dyn. 

The most likely cause of the force that drives the convectional motion is believed 
to be the buoyancy force which is given as pagéT, where a is the coefficient of cubical 
expansion, g is the acceleration of gravity and dT is the temperature departure from 
the surroundings. Putting a=4.5x10-* °C"! and g=800 cm/sec’, the buoyancy force 
becomes 0.0407 as has been shown by Bullard and Gellman“), For the geomagnetic 
dynamo, it has been also shown by them that there is evidence to suppose that dT 
is of the order of 10->~ 107°C, so that the buoyancy force would be of order 107° ~ 
10-*dyn. It would not be also impossible to suppose that dT is a little higher for 
localized eddies considered in this paper. 

Although the estimate above is very crude, we may expect that the pondermotive 
force produced by the interaction between the induced electric currents and the station- 
ary magnetic field could be of the same order of magnitude as that of the driving 
force. Since the pondermotive force is directed so as to counteract the fluid motion, 
the growth of the induced magnetic fields would become prohibited sooner or later 
because the motion is necessarily slowed down after some time when the pondermotive 


force is going to overcome the driving force. Whether such a process leads to an 


oscillation is not clear until a more detailed account of magneto-hydrodynamics becomes 
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available. The fact that magnetic fields of many types are induced as we have seen 
in Section 2 makes the situation more difficult. 

In the light of the above discussion, it seems doubtful that we have a monotonous 
approach to the steady state as shown in Fig. 4 in actual cases though the basic idea 
of induction model would still work in accounting for geomagnetic non-dipole fields 
and their secular variations. 
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Abstract 


The distribution of the sound velocity is obtained by the distribution of the 
absolute temperature and wind. When the direction of the sound is toward 
west at the middle latitude, the minimum velocities in summer and winter 
are estimated to be 210m/s at 70km in height and 280m/s at an altitude of 
about 90km respectively. The typical paths of sound waves will be shown by 
Snell’s law in the channel with the axis at these altitudes. The maximum of 
the sound intensity at an altitude of about 90km is reduced to about 0.055 dyne/cm*. 
The variation quantity of the density which is originated from sound is less than 
the value of density under the state of equilibrium. When the value of the 
density-variation reaches about 10-§ gm/cm® the wave form of the sound will be 
transformed into the wind at the altitude of about 90km. During the winter, at 
the middle latitude, a period of sound waves through the ionosphere comes to 


about 1 sec, corresponding to wave-length of about 280m. 


1. Introduction 


The phenomena of the sound propagation in the upper atmosphere have been 
investigated by E. Schrédinger, B. Gutenberg and others. In the upper atmosphere, it 
has been found that the velocity of sound has a minimum at an altitude of about 
80 kilometers in accordance with the distribution of temperature only. The propagation 
velocity of the sound wave is a function of the absolute temperature and affected by 
wind in the space. In this paper, we will consider the existence of the sound 
channel in the ionosphere and discuss the intensity and period of the sound in these 
channels. The main argument of this explanation is based on the results of the 
following articles, i.e. sound velocity, attenuation of the sound and wind speed in the 
ionosphere. 

The formula of the sound velocity in the upper atmosphere was obtained by M.J. 
Ference as T=3.488 10°C? (km sec). As to the attenuation of the sound in this 
altitude, E. Schrédinger introduced the terms of mean free path of the molecule, 
internal friction and heat conduction. According to the Rocket observation at the 


middle latitude, the west-component of wind speed at an altitude of 80km is about 
50 m/s in winter. 


( 54 ) 


The transmission of sound waves in the ionosphere is dependent on the vertical 
distribution of the sound velocity, which is detemined by the temperature and wind. 
The distributions of temperature and wind speed are given in Figs. (1) and (2). Based 
on the values of temperature and wind speed, the distribution of the sound velocity is 
shown in Fig. (3). It is assumed, here, that the direction of sound propagation is 
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2. Sound Transmission in the Ionosphere 
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Fig. 1. Temperature distributions in the upper atmosphere. 
(after W. Nordberg and W.G. Stroud) 
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\ Fig. 2. Wind velocity at the middle latitude. (from 
Hand Book of Geophys. United States Air Force 
Geophys. Research Directorate, Chapt. 5, p. 80) 
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Fig. 3. Distribution of the westward sound velocity in the upper atmosphere. The values 
of sound velocity is obtained by C=16.9 ./f77-+ (E-W component of wind velocity) 


westward and therefore the west component of wind speed adds to the sound velocity. 

We can then obtain the mode of the long propagation of the sound wave which depends 
upon the vertical distribution of the sound velocity. Snell’s law, ie. cos#=c/co is used 
to find ray paths of the infra-sound in the upper atmosphere. It is assumed that the 
ray path is circular and the velocity C is given by C=C, (l4rZ), where *tand is 


the radius of curvature of the ray, Z is vertical distance from the altitude of the 
minimum sound velocity and Cy) is the minimum sound velocity. Then 7 can be 
obtained as r= oye, We can get the path of the sound waves together with the 


value of radius and the critical angle which is obtained by the ratio of the minimum 
and the maximum sound velocities. 
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Fig. 4-a. Ray paths in the sound channel of ionospheric D-region. 
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Fig. 4-b. Ray paths through the E-region at the middle latitude in winter. The critical angle 0=33° 


and the horizontal distance of the propagation from the source to the apices are 30.72 km 
and 111.13 km. 


The typical paths of sound waves in the upper atmosphere thus obtained by the basic 
equation of ray acoustic are shown in Figs. 4-a,-b. 


3. Sound Intensity in the Ionosphere 


The energy of the source in producing sound waves is proportional to the product 
of acoustic impedance and the square of particle velocity. It is also proportional to 
the square of the pressure amplitude ». By numerous applications of the above two 
representations of sound energy and the relation of CAp/ pot between particles 
velocity and sound velocity, the relation of 4p and p are to be found. The intensity 
of the sound waves will be affected by the density of the medium. The sound waves 
with infra-frequencies will be transmitted through the ionosphere. The to-and-fro 
vibrations of the medium are transmitted through the layer with variable sound velocity. 
As shown in Fig. (5) the value of the density-variation 4p reaches about 0.055 dyne/cm”, 
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Fig. 5. Relation of the pressure amplitude P and the 
density-variation 4p of the sound waves. The maxi- 
mum intensity of the sound waves at the altitudes 
with various densities is shown. The curve of larger 
velocity is arranged on right side. 
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the wave form of the sound will be transformed into the current of the wind at the 
altitude of about 90km. It is possible when the value comes to 0.3 dyne/cm*, we can 
see the same phenomenon about 70km in height. If the pressure amplitude of 0.1 mb 
is detected on the ground, the sound wave will be not transmitted through the iono- 
sphere. The values of density py) at the altitude of the axis of sound channel in 
the upper atmosphere are about 10-* gm/cm’, 107? gm/cm* and 10-* gm/cm! respectively. 
When the sound occurs in the channels, it will be transmitted with lower density- 


variation than above mentioned densities. 


4. Period of the Sound waves in the Ionosphere 


The attenuation at the lower frequencies is smaller for the reason that the absorp- 
tion constant is proportional to square of frequency. The energy loss rates per km 
path length in an isothermal atmosphere at an altitude of 70km are 0.407 and 0.999 
at the wave-length of 14m and 200m respectively. The rate at 90km in altitude is 
shown in Fig. (6) under the condition of temperature of —45°C. Attenuation constant 
is given as 30.1//i® by E. Schrédinger, where / is mean free path of molecules and 4 
is a wave-length. From the relation of decreasing intensity with the wave-length, 
possible sound waves which propagate a long distance in the ionosphere is concluded 
to have larger periods than about 1 sec. If the source of the infra-sound is generated 
at the neighbourhood of the channel-axis due to the meteors, its part of long periods 
propagates into the sound channel of the ionosphere in wind direction. 


= 
S 


at about 90 &m in height 
---- a 70 km ” 


Energy rate per fm 
9 
Q 


(a) 100 200 300 
Wave length (m) 
Fig. 6. Relation of decreasing rate of the sound intensity 
with the wave-length at abont 90km in height. 
(after E. Schrodinger) 


5. Cloud of the Sound in the Ionosphere 


Attenuation of the sound must be small in order that the sound may exist in 
ionosphere. The propagation of the sound will be made easy by the actual existence 
of the sound channel. The density of the medium in the ionosphere is rare compared to 
that at low altitude, so the ordinary sound will disappear quickly. The sound that 
comes through the rare medium has long wave-length. We can find out the existence 
of the sound of 1 sec period from the relation of the wave-length of the sound with 
decreasing energy rate of the sound. The sound which has longer periods than 1 sec 
will exist easily at these high altitudes. The propagation distance of these sound 
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waves of the long-period is obtained by the calculation of an attenuation constant. 


The classical calculation of the absorption constant of the ordinary sound was indicated 


6.6 X . 
as aa ° Eo by Stokes, were C, 7, 9) and f are components of sound velocity, viscosity, 


density of the medium and frequency of the sound respectively. The formula of 
Schrodinger corresponds to that of Stokes mentioned above as to the plane wave of the 
sound. The mean free path of the molecules in connexion with the attenuation con- 
stant of the infra-sound in isothermal air was introduced by Schrodinger. The values 
of / at the standard temperature are about 0.075 cm and 0.7 cm at the altitudes of 70 km 
and about 90 km respectively according to the table of National Advisory Committee 
for Aeronautics. For example, if the wave-length are 210m and 280m at the 
altitudes of 70km and 90 km respectively, ie. the periods of the sound are 1 sec in 
both altitudes, the values of the attenuation constant under the standard temperature 
will be about 5.1x10-° and 2.7x10- at these altitudes. The decreasing rates of the 
sound-intensity, and not the intensity of shock waves, for distance of km are shown 
as the exponential function in Fig. 7. When the infra-sound takes the suitable path 
in the sound channel, it has the propagation of long distance. When the sources of 
the sound are created everywhere in the sound channel of E and D-region, the channel 
is filled with the cloud of the above-mentioned sound. In Fig. 7, the pressureampli- 
tude of the infra-sound at about 200km distance of travel in E-region is eotimated to 
be 76 % of its amplitude of near the source where it begins to be transformed into sound. 
The value of this pressure amplitude is about 0.042 dyne/cm?. 
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Fig. 7. Decreasing rates of the intensity of infra-sound as plane waves. 
On the other hand, Spherical waves spread in an infinite homogeneous 
medium in accordance with the law of inverse squares. 


6. Conclusion 


The possible axis of the sound channel in the ionosphere is taken to be the altitude 
of about 70km and 90km in summer and winter respectively. It is possible that the 
density waves in near ionosphere are inhaled into the sound channel. It can be expected 
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that the sound waves with pressure amplitude of about 0.055 dyne/cm? exist at the 
middle latitude of an altitude of about 90km. It may be considered that when the 
sounds of large period exists in above-mentioned channel, of some kinds of geomagnetic 
variation may occur according to these sound waves. Finally, the author wishes to ex- 
press his gratitude for the kind advice received from Honorary Prof. Dr. M. Hasegawa 
and Prof. Dr. M. Ota, University of Kyoto. 
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Letter to the Editors 


Effect of Solar Cosmic Rays on the Meteorite Radioactivity 


In recent papers” the spatial distribution of cosmic rays in the interplanetary 
space has been investigated through the production of radioisotopes in meteorites. 
Comparing the radioactivity of *7A and that of 389A, Stoenner, Schaeffer and Davis” 
concluded the spatial constancy of cosmic rays, whereas Fireman and DeFelice” re- 
marked a higher flux at distances closer to the sun. In this note it is pointed out 
that the contribution of solar particles to the production of radioisotopes which may 
favour the latter conclusion. 

The results of the above two measurements are not much different, as shown in 
Table 1, but their conclusions are opposite. 


Table 1. Ratio of 7A and 9A 


Bombarded by cosmic rays. Bombarded by protons 
Ref. 1. 2.0+0.3 1.5+0.2 (3 GeV) 
Ref. 2. 2.340.2 1.2+0.3 (2 GeV) 


The meteorite used for obtaining the above results is the same, the Hamlet chon- 
dritic meteorite which fell on October 13, 1959. About three months before this fall 
occurred the so-called July 1959 event which consisted of successive bursts of intensive 
solar protons. The integrated proton flux during this event is estimated to be as high 
as 4x10® protons/cm? for energies higher than 100 MeV.* This is higher than the 
Galactic cosmic ray flux during three months. j 

The amount of *“A of the mean lifetime t=49 days produced by ee irradiation 
by a constant flux J is proportional to r/, while that by an additional flux 4/7 during 
time T having occurred ¢ days ago is proportional to 4/-T exp (—t/r). Since J is 
2z times the unidirectional intensity, 0.18 cm~? sec™ sterad“!,” we have Jr~5x10"/cm?. 
The intensity of the solar particles at the meteorite, which was located far from the 
earth at the time of the July event, may be much smaller than that at the earth. 
Moreover, the steepness of their energy spectrum has to be taken into consideration. 
It may not be far from reality, on account of the above considerations, if we assume 
the effective flux of the solar particles responsible for the production of *’A in the 
meteorite to be about one tenth of the flux at the earth, 4/-T~4x10'/cm*. Thus we 


obtain 
A1-T exp(—t/t)/It=7 X 10°/5 X 107=0.1. 
This estimate does not claim that the *7A/°A ratio of the Hamlet meteorite should 
be 10% higher than that expected from the constant irradiation, but merely indicates 


that the contribution of solar particles is so appreciable that the observed results could 
be interpreted in this way. We would further like to point out that the meteorite 
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activity gives some information on the propagation of solar particles in a region far 
from the earth and suggests an experiment for their detection by means of artificial 
space probes. 


By Satio Hayakawa 
Physical Institute, Nagoya University. 
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